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INTRODUCTION 


This  report  summarizes  work  to  cate  on  a  new  method, 
simple  yet  accura:e,  for  analyzing  the  turbulent  boundary 
layer  for  the  variation  of  surface  shear  stress  and  heat 
transfer  under  fairly  arbitrary  conditions.  Earlier  work 
was  concerned  with  two-dimensional  compressible  flow  and 
was  summarized  in  the  report  by  White  end  Christoph  (1970) 
and  a  subsequent  published  paper.  White  and  Christoph 
(1972).  Another,  more  comprehensive  paper,  including  such 
effects  as  roughness  end  wall  transpiration,  has  also  been 
published:  White  and  Christoph  (1971).  The  latter  paper 
is  also  confined  to  two-dimensional  flow. 

The  present  report  extends  this  new  theory  to  ♦wo  new 
types  of  conditions:  1)  a  skin  friction  theory  for 
axisymnsetric  and  three-dimensional  flow  conditions;  and 
2)  a  new  analysis  for  two-dimensional  surface  heat 
transfer-  The  approach  retains  the  earlier  philosophy  of 
using  only  wall-related  temperature  and  velocity  profiles 
to  develop  integral  relations  for  skin  friction  and  heat 
transfer  which  are  devoid  of  shape  factors  or  integral 
thicknesses.  Thus,  for  example,,  the  axisymmetric  theory 
results  3n  a  single  differential  equation  for  wall  skin 
friction,  with  no  other  variables  appearing. 

It  is  convenient  to  di.ic®  the  report  into  tnree 
chapters.  The  first  chapter  develops  a  two-dimensional  thee... 
for  surface  heat  transfer  (or  Stanton  number)  for  compress¬ 
ible  flow  with  arbitrary  wall  temperature .  This  extends  the 
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earlier  report,  White  and  Christoph  (lS^O),  which  consider¬ 
ed  -kin  friction  only  and  eliminated  the  temperature  through 
a  Crocco  approximation.  The  new  theory  results  in  two  first- 
order  coupled  differential  equations  for  skin  friction 
coefficient  and  Stanton  number.  The  theory  is  compared  with 
variable  wall  temperature  end  pressure  gradient  experiments, 
and  an  appjoxim&te  correlation  is  found  which  relates  the 
Reynolds  analogy  factor  (c^/c^)  to  the  local  pressure 
gradient. 

The  second  chapter  develops  an  analysis  for  compress¬ 
ible  axi syrae t ri c  flow  with  pressure  gradients.  The  boundary 
layer  is  assumed  to  be  thick,  so  that  transverse  curvature 
is  important,  8nd  temperature  is  eliminated  through  the 
Crocco  approximation.  The  result  is  a  single  first-order 
differential  equatior  for  the  skin  friction  coefficient. 
Comparison  is  made  with  experiments  for  supersonic  flow 
along  thin  cylinders  and  bodies  of  revolution.  The  transverse 
curvature  effect  can  be  quite  large,  causing  skin  fraction 
increases  of  10O^  or  more  in  both  subsonic  and  supersonic 
flows.  Some  earlier  work  on  this  type  of  analysis  was 
given  in  a  paper  by  White  (1972). 

The  third  and  final  chapter  is  concerned  with  three- 
dimensional  incompressible  turbulent  skin  friction.  The 
distinguishing  characteristic  of  a  th  r  e  e  -  d  imen  s  ionc\l 
boundary  layer  is  a  crossflow  normal  to  the  freestream 
direction,  resulting  in  a  vector  wall  shear  stress  which  has 
both  streamwise  and  crossflow  components.  Since  both  shear 
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components  vary  with  both  surface  coordinates,  partial 
differential  equations  are  inevi table.  The  present  analysis 
uses  both  a  streaawise  ana  a  crossflow  1 aw- o  f - tb  e  ~ wall  to 
derive  two  first  order  coupled  partial  differential 
equations  for  the  two  surface  shear  components.  As  with 
previous  work,  there  are  no  shape  factors  o r  integral  thick¬ 
ness  present.  Also,  the  partial  differ''". Ma-*  equations 
are  of  such  simple  f o  m  ’hav  tht V  .  °f  xl'r  to  e 
simple  finite  difference  technique.  l’><  -  -tiers  :-."i 

applied  to  several  experimental  flows.,  am',  the  agreement  is 
very  premising,  particularly  considering  the  r  lative 
inaccuracy  of  present  three-dimensional  flow  «jasurements. 
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lhaoter  One 


Heat  Trans fex  Analysis 


1.1  Introduction 

The  purpose  of  this  chapter  is  to  develop  a  method  for 
calculating  heat  transfer  coefficients  for  arbitrarily  high 
Mach  numbers,  variable  wall  temperatures,  and  variable 
streamwise  pressure  gradients.  The  present  analysis  is 
restricted  to  steady,  two-dimensional  flow  of  a  perfect  gas 
in  a  compressible  turbulent  boundary  layer.  Hie  lateral 
pressure  gradient  will  be  neglected  since,  although  the 
lateral  gradient  strongly  influences  the  momentum  thickness 
in  high  Mach  number  flews,  it  appears  to  have  negligible 
effect  on  the  approach  proposed  here.  Two  coupled,  first 
order,  ordinary  differential  equations  are  derived  which 
have  for  dependent  variables  only  the  skin  friction  coeffi¬ 
cient  and  the  Stanton  number.  Stanton  numbers  are  calculat¬ 
ed  and  compared  to  the  data  of  Moretti  and  Kays  (1965)*  A 
relationship  among  skin  friction,  heat  transfer,  and  pressure 
gradient  is  presented. 

1.2  Development  of  the  New  Method 

The  equations  governing  compressible,  two-dimensional 
turbulent  boundary  layer  mean  flow  analyses  were  apparently 
first  given  by  Young  (1953)  In  the  following  form: 


4 


tL  d 


There  are  six  unknowns  (p  ,u, v,ho,q,r)  and  only  four 
equations,  so  that  further  relations  are  needed.  The 
finite  difference  methods  model  the  eddy  viscosity  and  eddy 
conductivity,  thus  correlating  the  variables  q  and  t  with 


local  conditions,  and  then  attack  the  complete  governing 
equations.  The  so-called  von  Karman  integral  methods 
eliminate  p  v  f ram  equations  (1-2)  and  (1-3)  by  equation 
(1-1)  and  then  integrate  with  respect  to  y  across  the  entire 
boundary  layer,  thus  obtaining  the  von  Karman  integral 
relations.  Additional  relations  for  these  integral  methods 
are  then  brought  in  as  correlations  between  integral  param¬ 
eters  . 

The  compressible  flow  analysis  of  White  and  Christoph 
(1970,  1972)  was  based  on  Crocco's  approximation  for  the 
energy  equation  -  cf.  Schlichtlng  (I960): 

p 

Prandtl  number  =  unity;  T  =  a  -*•  bu  +  cu  ,  (1-6) 

where  a,b,c  are  determined  from  boundary  conditions. 
Basically,  the  analysis  of  White  and  Christoph  (1970.1972) 
non-dimensionalizes  the  velocity  u  with  respect  to  wall 
quantities  and  then  integrates  equations  (1-1)  and  (1-2) 
with  respect  to  the  law-of-the-wall  variable  y+  =  yn*/  \>v, 
where  u*  =  (-r^/p^)^2  is  the  friction  velocity.  With 
Crocco‘s  approximation  one  can  compute  local  heat  transfer 
only  through  a  Reynolds  analogy.  Crocco’s  approximation 
will  not  be  made  here.  Instead,  the  energy  equation,  as 
given  by  equation  (1-3),  will  be  used.  The  velocity  u 
and  the  temperature  T  are  non-dimensionalized  as  follows: 


(1-7) 


uf  =  u(x,y)/u*(x)  and  T~  •=  T(x,y)/Tw(x)  . 


Quantitative  expressions  for  u"r  and  T+  are  needed. 
The  siinpie  incompressible  formula 


ln(y‘ ) 


O.  ~  G 


(1-8) 


and  a  similar  expression  for  Tf  are  not  adequate,  oecause 
we  need  to  include  compressibility,  heat  transfer,  and 
pressure  gradient  effects  in  th'  velocity  and  temperature 
correlations.  An  analytical  approach  of  finding  expressions 
for  u+  and  T+  was  desired,  so  that  as  little  empiricism  as 
possible  would  be  introduced  into  the  present  method.  The 
scheme  adopted  was  based  on  the  eddy  viscosity  approach  of 
Deissier  (1959) *  which  assumed  that  the  Prandtl  mixing 
length  approximation  could  b€:  extended  to  the  variable 
density  case  with  no  further  changes.  Let  us  first  derive 
an  expression  for  u  .  'Ihe  total  shear  stress  is  related 
to  an  eddy  viscosity  and  mixing  length  as  follows: 


T 


> 


(1-9) 


where  e  =  p«c2y^i|^i 


;  and  x  =  C.4  is  von  Karman's  constant.  Also,  near  the  wall, 

\ 

i 
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the  boundary  layer  is  approximately  a  Couette  flow  with 
negligible  acceleration,  so  that  equation  (1-2)  becomes 


(1-10) 


Non-dimensionaliring  equations  (1-9)  and  (1-10)  with  respect 
to  wall  variables,  we  neglect  the  laminar  portion  of  the 
total  shear  in  equation  (1-9)  (This  step  greatly  simplifies 
the  final  form  of  u^  but  has  a  negligible  effect  on  the 
numerical  values  of  u  . ) .  Equating  equations  (1-9)  and 
(1-10)  results  in 


=  (t*(i  ±  a  y*))1/2  (1-11) 

dy+  <  y+ 


where 


a 


v  , 
_ *  d£ 

Tu*  lx 
w 


(1-12) 


This  is  the  same  parameter  a  ,  termed  the  "pressure 
gradient"  parameter,  which  appeared  in  White's  (1969) 
analysis,  and  which  was  suggested  by  the  work  of  Mellor 
(1966).  It  is  an  ideal  parameter  for  this  method  because 
it  is  directly  related  to  the  shear  stress  without  any 
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integral  thicknesses,  we  now  need  an  expression  for  T  . 
Following  the  same  approach,  the  total  heat  flux  is  related 
to  an  eddy  conductivity  as  follows: 


p  /  x  ,  dT 
q  -  tt;  (  M-  +  e) 


(1-13 


where  F^.  is  the  turbulent  Prandtl  nu-.er.  Near  the  wall, 
one  obtains  from  equation  (1-3) 


q*  =  ut  +  q 


(1-14, 


Again,  non-dimensionalizing  equations  (1-13)  and  (l-l4), 
neglecting  the  laminar  portion  of  the  total  heat  flux  in 
e auction  (1-13),  arid  equating  equations  (1-13'>  and  (1-14) 
results  in 


o-gyu+d  +  g  y+))T+  1/2 

*y+  (1  +  o  y+)1/i? 


(1-15) 


where 


e  - 


PT% 


u*cT 
w  p  v» 


Y  =  ?rr 

P  w 


(l-lc> 


The  parameters  £  and  y  will  be  referred  to  as  the  "heat 
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transfer”  parameter  and  the  "compressibility"  parameter, 
respectively.  The  effects  of  the  parameters  a  ,  3  , 

and  \  on  u+  are  shown  in  Figure  (1-1).  Positive  a 
(adverse  pressure  gradient)  and  positive  3  (cold  wall 
heat  transfer)  raise  a+  above  the  incompressible  logarith¬ 
mic  law,  and  negative  a  (favorable  pressure  gradient) 
and  negative  3  (hot  wall  heat  transfer)  have  the  opposite 
effect.  The  parameter  y  ,  which  is  always  positive, 
lowers  u+  ceiow  the  incompressible  result.  If  more  than 
one  of  these  parameters  are  nonzero,  the  effect  is  roughly 
an  additive  one.  Precisely  these  affects  have  been  found 
experimentally,  e.g.,  Kepler  and  O'Brien  (1959),  Lee  et  al 
(1969),  and  Brott  el  al  (1969).  The  "heat  transfer” and 
"compressibility*'  parameters  have  the  same  qualitative 
effect  on  as  they  do  on  u+.  The  "pressure  gradient” 

4. 

parameter  alone  does  not  affect  T  .  However,  a  positive 
a  combined  with  either  a  positive  3  or  a  negative 
3  depresses  the  3  effect  on  T+.  A  negative  a  has 
the  opposite  effect.  The  compressibility  effect  is 
increased  with  a  positive  a  and  decreased  with  a  nega¬ 
tive  3  .  These  results  are  shown  in  Figure  (1-2). 

Equations  (1-11)  and  (1-15)  give,  for  zero  pressure 
gradient , 

=  =  3-2  v  u+  (1-17) 

du+  du+ 
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Figure  1.1.  ILLUSTRATION  OF  THE  EFFECT  OF  VARIOUS 
PARAMETERS  ON  THE  VELOCITY  LAW-0 F-THE-WALL 


FOR  COMPRESSIBLE  TURBULENT  FLOW. 


log  (y  'i 


Figure  1.2.  ILLUSTRATION  OF  THE  EFFECT  Of  VARIOUS 
PARAMETERS  ON  THE  TEMPERATURE  LAW-OF-THE-WALL 
FOR  COMPRESSIBLE  TURBULENT  FLOW. 


I 
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Integrating  equation  (2-17)  fives 


T  =  1  + 


0  u+  - 


v  u 


(1-16) 


which  is  the  form  of  Croceo’s  relation  used  by  White  and 
Christoph  (1970,  1972).  If  a  does  not  equal  zero,  one 
then  obtains 


ST*  _  ft  -  2  y  u*(l  ±  a  y*) 
Su*  1  f  a  y+ 


(1-19) 


which  is  not  a  Crocco  relation. 

Let  us  continue  the  method  of  solution,  by  first 
defining  a  compressible  stream  function  t  as  follows 


W 


P  u; 


hi 

Sx 


P  v 


(1-20) 


Eliminating  p  v  from  equations  (1-2)  and  (1-3)  by 
equation  (1-20),  and  non-dimensi<  nalizing  u  and  T  as  in 
equation  (1-7)  gives 


p  u*u+  |^(^+) 


ax 


(u*u+)  = 


u* 

V 

w 


Sr 

Sy+ 


(1-21) 
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and 


(  p  u*u+  5 - -I*-  ~r)  (c  T  T+  +  u-*^~)  = 

v  H  ox  ax  v  >  +'  v  P  w  2  ; 

w  oy 


u*  a _ 

V  -V  + 

v  ay 


(ut  +  q) 


(1-22) 


The  pressure  gradient  in  equation  (1-21)  is  related  to  the 
freestream  conditions  by  the  inviscid  Bernoulli  relation: 


dUe 

dx 


(1-23) 


The  y+-derivatives  in  equations  (1-?1)  and  (1-22^  are  no 
problem  because  these  equations  are  going  to  be  integrated 
with  respect  to  y4",  but  the  x-derivatives  must  be  calculated. 
The  functional  relationships  of  u^,  T+,  and  t  /i  (  ♦  has 
dimensions  of  viscosity)  are 

u+,  T+,  */(iw  -  fcn(y+,  a  ,  3,  y)  (1-24) 

The  r-derivatives  are  calculated  by  using  The  chain  rule: 


+ 


&  a 

35T  3o" 


+ 


a 

35T  oP 


.  dv  a 
+  “5x“  Ty- 
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For  skin  friction  ar.a  heat  transfer  calculations,  it  is 
convenient  to  define  further  non-dimensionalizations : 


x*  =  x/L  ;  7  =  Ue/U0 


\1  =  (2/Cf) 


1/2 


(1-26) 


=  T  J  T 
vK  wo 


X2  =  ^wTe/(pwU;Tw) 


where  L,  UQ,  and  are  constant  reference  values.  Now 
use  equations  (1-26)  in  the  definitions  for  y+,  a,  0,  and 
y  and  carry  out  the  partial  derivatives  in  equation  (1-25)' 
In  so  doing,  derivatives  of  v  ,  p..,  and  T  will  arise. 
Eliminate  such  derivatives  in  favor  of  (dT^/dx*)  and 
(d7/dx*)  through  equation  (1-J)  and  a  viscosity  power-law 
assumption  h*Tt.  where  x  is  taken  here  as  0.67  for  air. 
Also  use  the  relation  between  velocity,  pressure,  and 
temperature  in  the  freestream: 


e  v 


IT  =  f-^  '  ~  cp/cv)  *  (2-27) 

e 


The  net  result  is  that  equation  (1-25)  can  now  be  written  as: 
X5V  ^  dx*  Ty*  +  3:1  h  +  3  i  -  2v 


+  TV  V 

V  y  ~y 


_  X  x 

23  -rjr  +  2v 


,  X  3  (1/V)" 

>v.  *1  FL  *2 


(1-28) 


•f  -y*  (  j  +  i/2)  Tv4-  +  a(r'  +  i/2)  la  -  %  |gr_ 


1  14  6  V  +  <*£>*?  yr«*  4-  -  v4-  - 

X,  dx*  _  *0  2  e  7  **  '  ey  dp  Y 


*V- 


-* 
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where  the  primes  denote  differentiation  with  respect  to  x*. 
Equation  (1-28)  is  now  substituted  into  equations  (1-21) 
and  (1-22)  and  the  resulting  equations  integrated  from 
v '  =  0  (t  =  -  \  to  y+  =6  (t  =  0).  The  result  is 

W 


dx*'  (^-^Hj) 


(j.//)  H0  X -j  dx.^ 


X2  *3 


+  X1  7  'Xl‘6  'Hl"2H3'(iTi)  K  (»VVH^ 


(1-29) 


X1  t  3('r+1)  :VH4_  =  ‘V 


X14(1/V),,K6  dx2 


x; 


X1  7  V?E7+(^)Me  (3K6-h5-H^ 


(1-^0) 


-  X^  a(fr+l/2)H^+Hg  —  JL^('e/*w)  X2X1~"V^ 


where  is  an  effective  Reynolds  number  defined  as: 


UoL  ue  Te  1/2 

\ =  -r-  «r>  (’r) 

e  w  w 


(1-31) 


The  coefficients  are  integrals  involving  u^,  T+  and  the 
law-of-the-wall  parameters;  they  are  lifted  in  Appendix  A. 

It  is  possible  to  numerically  evaluate  these  H's  once 

4*  4*  4-  4- 

^  ^  >u  bn  >u  *T 

expressions  for  the  derivatives 
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~jr  ,  and  are  found.  Expressions  for  th=>se  derivatives 
38  3v 

are  obtained  by  differentiating,  with  respect  to  a,  3,  and 

+  4- 

y,  the  previously  derived  integrals  for  u  ,  T  ,  and 
Eqns. (1-11, 15, 20)  respectively,  at  a  fixed  y  .  Integral 
expressions  for  these  partial  derivatives  are.  given  in 

-f-  +  -f  4. 

Appendix  A;  the  change  of  variables  du  =  (>u  />y  )dy  has 

been  made.  The  coefficients  are  then  evaluated  by 

integrating  the  relations  in  Appendix  A  from  u+  =  0  to 

u+  =  u+  =  \t(T  /T  )^2for  various  values  of  a,  3,  and 
e  *lv  e  w' 

* 

Y.  The  turbulent  Prandtl  number  is  taken  as  0-9  • 

x 

The  starting  value  of  T+  is  1.0,  and  the  initial  value  of 
y+  is  taken  to  be  0.1106  =  exp (-5 *5/2. 5)  to  match  with  the 
incompressible  logarithmic  law  (Eqn.1-8)  in  the  limit.  Curve 
fit  expressions  for  the  coefficients  are  also  given  in 

Appendix  A. 

Tiie  parameters  a,  3,  and  y  are  related  to  the  dependent 
variables  X^  and  as  follows: 

a  -  X3  (l/V)'/*^ 

YUe?  "  PT  Me  (Te/Tw)  (1"32) 

K  =  PT  X1  X2  >4  (H)(Te/Tw) 

* 

i 

* 

PT  =  0.9  was  found  experimentally  by  Simpson,  Whitten  and 
Moffat  (1970).  It  is  also  the  value  assumed  in  most 
finite-difference  calculations,  e.g.  Cebeci  0-9^1 )• 


17 


,1 


Equations  (1-29)  and  (1-30)  are,  then,  two  coupled, 
first  order,  nonlinear,  ordinary  differential  equations  in 
and  \ p  only-  The  freestream  Mach  number  and  wall  tempera¬ 
ture  distributions  must  be  known,  plus  the  values  of  X^  and 
X2  at  some  initial  x*.  Unfortunately,  5+  appears  in 
equation  (1-29)  and  must  be  calculated.  However,  6+  need 
not  be  Known  precisely  for  accurate  calculations  of  x^  and 
Xg.  An  approximate  way  of  calculating  6  ,  and  the  one 
recommended,  is  to  first  assume  a  Crocco  law  of  the  form  in 
equation  (1-10)-  With  this  expression  for  T+,  equation 
(1-11)  may  be  integrated  analy tic ally  as 


Pin*1  =  Sin*1  (^"o'6)  + 

^  FI 


+  ~i^2  (p-po>  +  1x1  ^7+TF^:) 

o 


(1-33) 


where 


q  =  (&2+ky)1/'2  and  P  =  (l+ay+)1//2 


Equation  (1-33)  needs  to  be  initialized.  By  plotting  u+ 

versus  y+,  as  in  Figure  1-1,  it  is  seen  that  the  curves  for 

various  values  of  a,  0  and  Y  converge  at  y+  slightly  less 

than  ten.  The  initial  conditions  u+  =  0  and  y+  =  0.1108 

o  o 

were  taken.  Unfortunately,  y  is  an  implicit  function  of 
u+  so  that  an  iteration  scheme  must  be  used. 

The  dependent  variable  x^  Is  directly  related  to  the 
Stanton  number 

X  (k-1)  Mf  (Te/V, 

ch  *  — v^r^r 
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where  the  Stanton  number  has  been  defined  as 


= 


% 


c  U  Tf  “) 
r  o  ev  aw  w' 


(1-35) 


Since  \1  is  defined  as  (2/c^.)  '  * ,  equations  (1-29)  and 
(1-30)  give  and  c.  directly  -  no  auxiliary  information 
is  needed  about  the  behavior  of  integral  thickness  or  shape 
factors.  Another  important  and  unique  property  of  equations 
(1-29)  and  (1-30)  is  an  explicit  flow  separation  criterion. 
As  separation  is  approached,  d\-j/dx*  approaches  infinity 
and  hence  c,.  approaches  zero  -  which  is  precisely  the 

x 

definition  of  two-dimensional  separation.  One  need  not 
worry  about  attempting  to  predict  separation  according  to 
a  particular  value  of  the  shape  factor. 

One  basic  assumption  is  made  in  the  development  of  the 
present  method;  that  is,  that  the  velocity  and  the  tempera¬ 
ture  profiles,  correlated  by  inner  variables  only,  are  valid 
all  the  way  to  the  edge  of  the  boundary  layer.  This 
assumption  leads  to  a  slight  miscalculation  in  6+,  as  seen 
by  Figure  1-3*  The  outer  wake  is  not  predicted  accurately. 

But,  as  was  pointed  out  earlier,  it  is  not  terribly 

+ 

imrortant  to  know  5  precisely  in  order  to  obtain  accurate 
skin  friction  and  heat  transfer  coefficients.  Also,  in  a 
supersonic  turbulent  boundary  layer,  the  wake  almost 
disappears.  Of  course,  there  will  be  some  error  in  the 
momentum  and  the  displacement  thicknesses,  if  one  wishes  to 
calculate  them.  In  passing,  it  should  be  noted  that,  from 
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- f. 

log (y  ) 


Figure  1-3.  COMPARISON  OF  THE  PRESENT  NEAR-WALL 


THEORY  WITH  THE  ACTUAL  LAW-OF-THE-WAKE . 


equation  {1-11 ),  if  ay*  is  greater  than  minus  one,  then 
^  ^  + 

•~f  !  —f-  |  <r  0.  Thus  the  velocity  gradient  is  negative,  or 
tne  velocity  insice  the  boundary  -3yer  has  become  greater 
than  the  freestream  velocity.  to  avoid  this  difficulty, 
it  is  suggested  that  one  take  6*  =  -1/a  in  the  case  of 

a  favorable  pressure  gradient.  This  is  not  believed  to 
handicap  the  present  methoc. 

1-3  Comparison  With  Heat  Transfer  Experiments 
1.3.1  Flat  Plate  Fata 

First,  it  seems  appropriate  to  consider  a  special  case, 
d?  dq 

VI  r  T) 

namely  when  ^  =  G.  For  this  case,  the 

present  analysis  reduces  to  the  analysis  given  by  White  and 
Christoph  (19FC,  19 Equation  (1-29)  reduces  to: 

d).-, 

Hi  37*  =  \  »  (1-36) 

which  integrates  to: 

0.455 

~  ,  — - -  r - - - n  a.,v- 

f  A- In*  0.06  (?7Twr^/A  "  J' 

where  A  is  the  van  Lriest  flat  plate  parameter,  given  by 
equation  (2-32)  and  plotted  in  Figure  2-3- 


I 


21 


litis  result  is  expressed  in  terns  of  a  flat  plate 
compressibility  transformation  for  skin  friction  as  follov/s 
(in  the  .notation  of  Spalding  3nd  Chi  (1964)): 


Cf  “  FcCflnc.<RxFib'> 


where  F  =  A 
c 


(1-36) 


1+7 


and  Pftt  -  (Te/Tw)-/A 


The  present  theory  and  five  other  theories,  the  five  most 
popular  and  accurate  kno;«n,  were  compared  to  427  adiabatic 
and  230  cold  wall  heat  transfer  data  points.  The  five 
other  theories  were  cast  in  the  same  form  as  equation  (I-38), 
but  of  course  and  were  different  for  each  theory- 
All  theories  were  computed  with  the  incompressible  skin 
friction  formula  of  Spalding  and  Chi  (1964),  which  apparent¬ 
ly  gives  the  best  agreement  with  incompressible  friction 
data.  Both  the  root  mean  square  error  and  the  mean  abso¬ 
lute  error  were  computed.  The  results  are  shown  in  Table 
1-1.  For  adiabatic  wall,  the  methods  of  Spalding  and  Chi 
(1964),  van  Driest  (1956b)  and  the  present  theory  are 
equally  accurate.  For  cold  wall  heat  transfer,  the  present 
theory  has  the  lowest  percent  error,  with  Moore  (1962) 
second.  The  present  method  is  thus  the  most  accurate  flat 
plate  theory  available  at  the  presen:  time,  with  the  added 
advantage  that  it  can  be  extended  to  more  general  freestream 
and  wall  conditions  (White  and  Christoph  1971). 
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TABLE  1-1 


COMPARISON  OP  SIX  THEORIES  WITH  PLAT  PLATE  FRICTION  DATA 


.ADIABATIC :  427  POINTS 

COLD  WALL:  230  POINTS 

AUTHOR 

RMS  <  ERROR 

ABS  %  ERROR 

RMS  %  ERROR 

- 1 

ABS  %  ERROR 

Eckert  (1955) 

12.44 

9.06 

_ 

29.49 

25.56 

Moore  (1962) 

8.67 

B 

17.69 

13.08 

Sommer  and 
Short  (1955) 

9.40 

7.77 

23.55 

20.14 

Spalding  & 

Chi  (1964) 

7.59 

5.46 

21.13 

16.94 

Van  Driest 
#2  (1956b) 

7.55 

5.45 

17-49 

13.81 

Present 

Theory, 

Eqn.  (1-37) 

7.80 

5-26 

14.31 

11.28 

3  -  1/2 

Note:  RMS  Error  =  (^-  I  e^)  ; 

MEAN  ABS  Error  =  ~  I  |  e± | 


Bata  Source:  White  and  Christoph  (1970). 


I 
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1.3.2  Data  of  Moretti  and  Kays 

Next,  the  present  method  was  compared  to  the  data  of 
Moretti  and  Kays  (1965).  Their  data  were  taken  in  a 
rectangular  duct  which  followed  a  system  of  rectangular 
nozzles.  This  test  section  was  designed  to  provide  various 
desired  distributions  of  stream  velocity  and  surface  tempera¬ 
ture.  The  velocities  and  heat  transfer  rates  for  these 
experiments  were  low,  so  the  boundary  layers  were  essentially 
constant  property  layers.  However,  the  important  effects 
of  wall  temperature  variation  and  f reestream  velocity 
variation  were  present.  In  fact,  predicting  heat  transfer 
under  varying  f reestream  conditions,  especially  accelerat¬ 
ing  flows,  has  become  of  practical  concern  in  such  problems 
as  cooling  of  gas  turbine  blades  and  rocket  nozzles. 

Figure  1-4  shows  the  Stanton  number  as  a  function  of 
position  along  the  rectangular  duct  for  a  constant  f ree¬ 
stream  velocity  and  for  a  sharp  decrease  in  the  wall  tempera¬ 
ture.  No  data  was  given  by  More  ,ti  ana  Kays  (19^5)  over  the 
first  one- third  of  the  flow  field.  Since  the  Stanton 
number  must  be  known  at  the  starting  point  of  integration  of 
equations  (1-29)  and  (1-30),  and  since  it  was  desired  to 
predict  the  effect  of  the  wall  temperature  decrease  which 
occurred  in  the  firs*,  one-third  of  the  flow  field,  equations 
(1-29)  and  (1-30)  we ' -e  integrated  forward  and  backward  from 
the  first  data  point.  Skin  friction  coefficients  were  not 
given  by  Moretti  and  Kays,  so  the  starting  value  of  was 
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chosen  such  "hat  c^  -  0.'--  c^.  The  curve-fit  wall  tempera¬ 
ture  distribution  was  chosen  to  be  of  the  fora 

m 

^  =  a  -  b  cos  tt(x*-x*)/c  (1-39) 

‘  wo 

where  (a,b, c)  were  fitted  constants  and  xj  was  the  x* 
position  where  the  vail  temperature  change  started-  The 
present  theory  and  the  digital  computer  finite  difference 
method  of  Herring  and  Mellor  (196c;  are  compared  to  the 
data  in  Figure  1-4.  Agreement  by  both  methods  is  excellent. 

The  second  Stanton  number  distribution  measured  by 
Moretti  and  Kays  consisted  of  a  wall  temperature  change  plus 
an  accelerating  freestream  velocity.  Again,  the  wall 
temperature  distribution  was  fitted  by  equation  (1-39) » 
but  the  freestream  velocity  distribution  was  fitted  by  a 
polynomial.  The  present,  method  and  the  method  of  Herring 
and  Mellor  are  compared  with  the  data  in  Figure  1-5-  Once 
again  both  methods  accurately  predict  the  data.  A  more 
severe  test  of  the  present  method,  and  of  the  method  of 
Herring  and  Mellor,  is  the  steep  favorable  pressure 
gradient  as  given  by  Moretti  and  Kays'  data  in  Figure  1-6. 

If  one  believes  in  a  constant  Reynolds  analogy  factor,  one 
would  expect  an  increase  in  the  Stanton  number  or  at  most  a 
constant  or  slightly  decreasing  Stanton  number.  As  seen  in 
Figure  1-6,  the  Stanton  number  decreases  considerably. 

Both  the  present  theory  and  the  method  of  Herring  and  Mellor 
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Figure  1.4.  COMPARISON  OF  THEORY  WITH  ZERO  PRESSURE 
GRADIENT  DATA  OF  JSORETTI  AND  KAYS  (1965) 
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Figure  1.7.  COMPARISON  OF  THEORY  WITH  ADVERSE  PRESSURE 
GRADIENT  DATA  OF  M0RETTI  &  KAYS  (1965). 
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qualitatively  predict  this  decrease.  However,  neither  the 
present  method  nor  the  method  of  Herring  and  Kellor  are 
able  to  fully  recover  after  the  favorable  pressure  gradient 
has  ended  -  the  present  method  doing  somewhat  better  than 
the  method  of  Herring  and  Mellor. 

The  effect  of  an  adverse  pressure  gradient  on  the 
Stanton  number  is  shown  in  Figure  1-7.  There  is  a  notice¬ 
able  difference  between  the  present  method  and  the  method 
of  Herring  and  Mellor  toward  the  end  of  the  flow  region. 

The  present  method  shows  a  slight  increase  in  the  Stanton 
number,  whereas  Herring  and  Mellor  show  a  slight  decrease. 
The  lata  of  Moretti  and  Kays  remain  approximately  constant 
and  are  in  reasonable  agreement  with  either  theory. 

l.A  Conclusion 

The  purpose  of  this  chapter  was  to  develop  a  simple 
and  accurate  method  for  analyzing  heat  transfer  and  skin  fric¬ 
tion  in  two-dimensional,  compressible  turbulent  boundary 
layers.  This  method  uses  law-of- the -wall  velocity  and 
temperature  correlations  which  include  pressure  graident, 
heat  transfer,  and  compressibility  effects.  After  integrat¬ 
ing  the  continuity,  momentum,  and  energy  equations  across 
the  boundary  layer,  two  coupled  differential  equations 
result  which  have  as  their  dependent  variables  the  skin 
friction  coefficient  and  the  Stanton  number.  Contrary  to 
the  von  Karmen  technique,  integral  thicknesses  and  shape 


factors  are  not  used.  The  present  heat  transfer  analysis, 
which  is  for  arbitrary  freestream  velocity  and  wall 
temperature  distributions,  compares  favorably  with  the  data 
of  Moretti  and  Kays  (1985)  and  is  as  accurate  as  the 
sophisticated  finite  difference  technique  of  Herring  and 
Mellor  (1968). 
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Chapter  Two 


Axisymmetric  Compressible  Flow 


2.1  Introduction 

In  this  chapter,  our  basic  approach  is  applied  to 
derive  a  method  for  computing  skin  friction  in  a  thick 
axisymmetric,  compressible  turbulent  boundary  layer.  The 
analysis  is  restricted  tc.  steady  flow  and  a  perfect  gas  is 
assumed  for  convenience.  The  freestream  Mach  number  and 
pressure  gradient  may  be  arbitrarily  large  and  variable. 
Nonadiabatic  wall  temperature  is  allowable,  but  the 
variations  should  be  "modest",  that  is,  the  streamwise 
derivative  of  wall  temperature  is  neglected.  The  tempera¬ 
ture  is  eliminated  in  favor  of  velocity  through  the  Crocco 
approximation.  It  has  already  been  shown  that  this 
simplification  has  little  effect  on  the  accuracy  of  skin 
friction  computation  (White  and  Christoph  1972). 

The  transverse  curvature  effect  is  introduced  through 
a  coordinate  change  suggested  by  Rao  (196? ) .  The  resulting 
differential  equation  for  skin  f  ration  (£q.  2-25)  Is  thus 
valid  for  arbitrarily  large  ratios  of  the  boundary  layer 
thickness  to  the  body  radius. 

Rao's  transformation  was  used  by  White  (1972)  to 
compute  the  skin  friction  in  turbulent  incompressible  flow 
past  a  long  cylinder  of  constant  radius.  The  basic  results 
of  that  paper  can  be  summarized  as  two  formulas  for  turbu- 
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lent  Incompressible  skin  friction  on  very  long  cylinders: 

r  "1 

C„  =  0.0015  +  0.20  +  0.016  (x/r  )0'4  R-1'3 

i  OX 

(2-1) 

=  0.0015  +  0.30  +  0.015  (L/ro)0*4  R£1'/3 

Here  rQ  is  the  cylinder  radius,  x  is  the  axial  distance 
from  the  leading  edge,  and  ?  =  Uqx/v  is  the  local  Reynolds 

number  for  an  assumed  constant  freestream  velocity  UQ.  The 
quantity  is  the  total  friction  drag  coefficient  on  a 
cylinder  of  length  L.  Equations  (2-1)  are  accurate  to  +  5# 
over  the  entire  turbulent  flow  range  (10^<Rj.  <  10^)  and  for 
cylinder  lengths  up  to  (L/r^)  =  10^.  Similar  results  will 
be  obtained  here  for  compressible  flow  along  a  cylinder. 

2.2  The  Axisymmetric  Law-of-the-Wall 

In  order  to  use  the  present  method,  it  is  necessary  to 
have  a  realistic  expression  for  the  law-of- the- wall  in  a 
thick  axisynmetric  boundary  layer.  This  law  differs 
considerably  from  the  two-dimensional  law  because  of  the 
effect  of  the  radial  variable  in  the  momentum  equation. 

With  reference  to  Figure  2-1,  let  the  body  radius  be  rQ(x). 
Within  the  sublayer,  both  the  convective  acceleration  and 
the  pressure  gradient  terms  in  the  momentum  equation  are 
negligible.  The  equation  reduces  to: 


or:  rt  =  r0  xw  =  constant  =  r  uw  (2-2) 

With  t  thus  approximated  by  the  viscous  stress  only,  we  may 
integrate  Eqn.  (2-2)  to  obtain  an  approximate  velocity 
profile  in  the  sublayer: 

“w  u  4  ro  TW  ln(r/r0)  (£-3) 

2 

If  we  now  define  the  shear  velocity  u*  such  that  T  =  p  u*  , 

W  n 

Eqn.  (2-3)  iaay  be  rewritten  in  law-of-the-wall  form: 

u+  =  r*  ln(r/ro)  ,  (SUBLAYER)  (2-4) 

Hh  -f- 

where  u  =  u/u*  and  r  =  p  u*r  /a  .  Equation  (2-4)  is  in 

u  w  u  w 

marked  contrast  to  the  two-dimensional  sublayer,  for  which 


u+  simply  equals  y+,  as  in  Chapter  1.  Since  y+  correlates 
the  entire  wall  region  in  two-dimensional  flow,  it  was 
suggested  by  Rao  (196?)  that  the  variable  r*  ln(r/ro)  in 
Eqn.  (2-4)  would  similarly  correlate  an  axisymmetric  flow. 
Thus  Rao's  hypothesis  for  incompressible  axisyrranetric  flow 
is  that 

u+  =  fcn(Y^),  where  Y  =  rQ  ln(r/rQ)  .  (2-5) 

Equation  (2-4)  would  hold  in  the  sublayer,  while  the  loga¬ 
rithmic  overlap  layer  would  be  characterized  by  the  relation 

u+  =  i  ln(Y+)  +  B  (  *=  0.40,  B  =  5-5)  (2-6) 

Rao  (1967)  and  later  rfhite  (1972)  showed  that  Eqn.  (2-6) 
is  an  excellent  approximation  for  wide  variety  of  thick 
axisymmetric  boundary  layers.  As  the  boundary  layer  becomes 
thicker,  this  axisymmetric  wall  law  becomes  valid  across 
the  entire  boundary  layer,  that  is,  the  outer  or  "wake" 
layer  becomes  vanishingly  small. 

Note  that  Eqn.  (2-6)  is  consistent  with  a  turbulent 
eddy  viscosity  e  and  mixing  length  l  given  by 

e  =  p  <  Yc  (r/ro)  ^  cos ( o )  ,  (2-7) 

or  i  =  < Y  (r/rQ)1/2  , 

where  o  =  tan-^(drQ/dx) .  For  this  report,  we  shall  take 
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d  =  0°,  that  is,  we  neglect  the  slope  of  the  body  and  take, 
approximately,  r  =  r  +  y.  The  error  is  negligible  except 
near  the  stagnation  point,  which  does  not  concern  us.  Since 
Y  as  defined  by  Eqn.  (2-5)  is  smaller  than  y,  t  as  defined 
by  Eqn.  (2-7)  is  smaller  than  its  two-dimensional  equivalent 
*  =  Ky»  the  physical  explanation  being  that  a  cylinder  has 

less  ability  to  create  turbulent  shear  than  a  plane  surface. 

Equation  (2-7)  may  be  used  to  derive  a  law-of- the- wall 
for  compressible  flow  with  pressure  gradient.  As  in  the 
two-dimensional  case  (White  and  Christoph  1571),  we  neglect 
the  convective  acceleration,  and  the  axisymmetric  streamwi^e 
momentum  equation  becomes 

7F<rT>  -  r^T  *  0  (2-8> 

Integrating  from  the  wall  (r=r  )  outward,  we  have: 


This  may  be  integrated  again  by  assuming  that  viscous  shear 
is  negligible  and  that  x  =  e  (vuAy),  where  e  is  given  by 
Eqn.  (2-7).  Substituting  for  x  and  rearranging  in  terms  of 
law-cf-the-wall  variables,  we  obtain: 


du 

dY< 


<Y 


(p 'l  +  f 


,  ov+/  'r  _l/2 

<  /r°- 1) 


where  a 


V  „  dPe 
7~v*  H3T 

W 


(2-10) 
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The  pressure  gradient  parameter  a  is  the  same  quantity 
which  appeared  in  previous  analyses  for  two-dimensional 
flow  (Chapter  1  and  White  1969). 

Equation  (2-10)  can  he  integrated  for  the  velocity 
profile  u+(Y+)  if  we  make  an  assumption  about  the  density 
variation.  Since  we  are  concerned  only  with  predicting 
skin  friction,  we  assume  a  perfect  gas, 

p  =  p  R  T,  or:  p/p  =  T/T  .  (2-11) 

2 

plus  a  Crocco  approximation,  T=a+bu+cu,  which  may 
be  rewritten  in  terms  of  wall  variables  as  follows: 

+  +2 

T/Tw  -  1  +  6  u  -  fu  =  p/p  ,  (2-12) 

where  0  and  y  are  the  same  parameters  used  in  previous 
work  (White  and  Christoph  1970,  1972): 

0  =  a  v/(Twkwtf*)  =  transfer  parameter; 

(2-13) 

A  p 

v  =  r  u»c/(  2c  T  )  =  compressibility  parameter. 

r 

Equations  (2-10)  and  (2-12)  may  now  be  combined  and  integrat- 
ed  to  obtain  the  desired  law-of-the-wall  u  (7  ,a,0,y).  Of 
particular  interest  is  a  closed  form  which  can  be  obtained 

for  the  special  case  (a=O,0=O):  j 

I 

u+(Y+,0,0,y)  =  -jL  sinr^  ln(Y+A*)  1  (2-14)  | 

where  we  take  Y^  =  0.1108  so  that  the  formula  reduces  in  the 

V 

limit  of  large  radius  to  the  two-dimensional  law-of-the-wall.  i 
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Equation  (2-14)  corresponds  to  supersonic  adiabatic  flow 
along  s  cylinder,  for  which  experimental  thick  boundary 
layer  data  is  available  (Richmond  1957). 

Some  velocity  profiles  obtained  from  integrating 
Eqn.  (2-10)  are  shown  in  Figure  2-2.  The  effect  is  generally 
the  same  as  in  the  two-dimensional  flow  study  (White  and 
Christoph  1970),  with  positive  a  (adverse  pressure  gradient) 
raising  the  curves  above  the  incompressible  log  law  and  y 
(the  compressibility  effect)  tending  to  lower  the  curves. 

Note  that  the  parameter  r*  has  no  effect  on  the  curves 
unless  a  is  finite. 


2.3  Derivation  of  the  Basic  Differential  Equation 

We  now  assume  that  the  compressible  law-of-the-wall 
for  a  thick  axisyuaietric  turbulent  boundary  layer  is  known 
as  the  integrated  combination  of  Equations  (2-10)  and  (2-12): 


.2  1 


u+(Y+,a,£,  X^q)  =  J  kYjt(1  +  0u+  -  Yu+  )2 
0 

Tl  -  aro  (l  -  e2Y+/ro)J  dY+ 


(2-15) 


~TT 


This  relation,  if  accepted  as  a  reasonable  approximation 
across  the  entire  boundary  layer,  provides  closure  to  the 
problem  of  computing  the  wall  shear  stress  on  a  body  of 
revolution.  That  is,  we  can  now  derive  a  single  different¬ 
ial  equation  for  Cj.(x)  under  arbitrary  flow  conditions. 

The  boundary  layer  equations  for  compressible,  axi- 
symmetric,  turbulent  boundary  layer  flow  are  given  by,  for 
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Y  =  r  ln(r/r  ) 
o  o 


SOME  EXAMPLES  OF  THE  AXI SYMMETRIC  LAV-OF-THE-WALL 
COMPUTED  FROM  EQUATIONS  (2-10)  and  (2-12). 


example.  Herring  and  Melior  (1968): 

a)  Continuity: 


+  |rr(Pvr) 


0 


(2-16) 


b)  Momentum: 


pur 


“(rT) 
>y v  ' 


pvr  — 


CD 

*  c 

dx" 


(2-17) 


c)  Energy: 


pur  *x~  + 


:»h 


*  r  ,  . 


Pvr  —  =  “y  r(<?  +  r)  (2-18) 


where  the  notation  is  the  same  as  in  Figure  2-1.  As  in 
Chapter  1,  the  quantity  hQ  *  h  +  u  /2  is  the  stagnation 
enthalpy.  Now  in  fact  Eqn.  (2-l6),  although  included  for 
completeness,  is  not  needed  here,  since  the  temperature  relates 
directly  to  velocity  through  the  Crocco  approximation, 

Eqn.  (2-12).  As  mentioned,  the  effect  of  this  simplification 
upon  skin  friction  is  slight,  but  knowledge  of  the  Jtanton 
number  distribution  is  lost.  If  desired,  one  could  extend 
the  temperature  wall  law  technique  of  Chapter  1  to  this  axi- 
synsnetric  case  for  computing  wall  heat  transfer. 

The  continuity  relation  (2-16)  is  satisfied  identically 
by  defining  an  axisyraaetric  stream  function  *  such  that: 


^  ■  0  u  r 


= 


pvr 


(2-19) 
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Since  p,  u,  r,  and  y  are  each  related  to  wall  variables,  it 
fellows  from  the  first  of  Eqns.  (2-19)  that  +  itself  is  a 
law-of- the- wall  variable: 

Y+  l 

j/uwro  =  j  £-.(!)*  ^  dY+  =  fcnrY+,a,P,Y,r^l  (2-20) 
0  Pw  ro 


This  fact  is  important  in  computing  the  velocity  v  from  the 
second  of  Eqns.  (2-19)-  Utilizing  wall  law  variables 
wherever  possible,  the  momentum  relation,  Eqn.  (2-17),  now 
becomes : 


,  -» 
pr  if  u 


=  -  r 


dPe 

35T 


(2-21) 


+  ^  Z—Avr 

vw  *V+V 


where  i#  =  is  the  wall  friction  velocity.  The  y+ 

derivatives  are  related  to  Rao's  variable  Y+  through  the 
relation 


(2-22) 


which  follows  from  Eqn.  (2-5)*  The  x  derivatives  must  be 
handled  by  the  chain  rule,  since  each  of  the  parameters 
(a,0,y,  r*)  in  the  law-of- tne-wall  is  a  function  of  x. 
Thus  we  substitute 


>X  ^y+ 


+ 


*£  2l_  +  2£iL_  +  hX 

7\x  -  >x  >v 


o 
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2 _ 


(2-2  3) 


In  carrying  out  the  partial  derivatives  in  Eqn.  (2-23),  we 
retain  our  previous  approximation  (White  and  Christoph  1970, 
7 972 )  of  neglecting  the  derivatives  with  respect  to  x  of 
pw,  T  ,  uw,  and  qw  wherever  they  appear  in  3,9,  y  ,  etc.  It 
is  felt  that  tnese  terms  have  only  a  small  effect  on  wall 
skin  friction. 

It  remains  only  to  carry  out  our  basic  integral 
procedure,  namely,  combining  Eqns.  (2-21,  22,  23)  and  inte¬ 
grating  the  entire  equation  with  respect  to  Y+  across  the 
boundary  layer  from  (Y^  =  0,  T  =  Tw)  to  (Y+  =  Y+,  r  =  0). 

The  result  is  a  single  first  order  differential  equation  for 
the  wall  friction  velocity  u*(x).  It  is  convenient  to 
define  dimensionless  variables  as  in  Chapter  1: 


X  =  >/27cJ  ;  x*  =  x/L  ;  V  =  Ue(x)/UQ  ,  (2-24) 


2 

where  C^>  =  2Tw/(peUe) ,  and  L  and  UQ  are  suitable  reference 
values.  In  terms  of  these  variables,  the  final  integral 
momentum  equation  becomes 


&>;*!  -  3°a2)  a,)  -  i  Ag(i/v)  =  aj_v 


-  A. 


1 

r 


4 

_  X 

o 


(2-25) 


1  1  ? 

where  ^  =  (UQL/ve)(uf  )(Te/Tj7  and  ?  =  ±  X 

r2(e?^-l). 
o v  7 


42 


Equation  (2-25)  is  the  central  result  of  this  chapter  -  a 

first  order  ordinary  differential  equation  for 

subject  only  to  some  assumed  initial  condition  Mx*0)  anc* 

Known  flow  conditions  U  (x).  T  (x),  T  (x),  r  (x),  and  M  (x). 

e  e  w  o'7  e '  7 

The  functions  (A^,  A?,  A-,)  arise  from  the  integration 
across  the  boundary  layer  and  are  defined  as  follows: 


A1 


}e{(Te/Tw)exp(2Y+/r^r  U+*  -  Bu+  |g+ +  2y  u+  |i+  ] 


_  3*  3u 

*wro  38  37+ 


2y  ) 

”wro  3Y+  f 


(2-26) 


k2~  !  (  e*P(2*’+/r*)  u+  -2—  — 


u  r 
w  o 


3*  3u+  \ 

1ST?-  r 


(2-27) 


In  order  for  Eqn.  (2-25)  to  be  viable,  it  is  necessary  to 
have  analytic  approximations  for  these  coefficient  functions. 
After  extensive  numerical  computation  of  Eqns.  (2-26,  27, 

28),  using  Eqn.  (2-15)  for  the  wall  law,  we  offer  the 
following  curve-fit  approximations: 


A1  " 


2  3 


A_  = 


c  ^ _ J  0.48  1+  0.2  S  +  0.0018  S£  \ 

8.0  expj  - nr'tTlZ  0T -  i  ’ 


(2-29) 


where 


exp(0.4A  /A)/r+  , 


and  where 


-  V</(*  l>  ’ 


jar*  (exp(2Y*/r*) 


^  _ /  0.84  x  +  0.2  s  +  0.0018  \  -,n\ 

0.066  expi  - A  '  n  -f  072  Ql -  V  (^-30^ 


=  0.1108  exp(  0.4 ^  -  2.2  Q2/3)  , 

where  S  =  sin’1  -f(2Yu+-  3)7(0^+  4y)1//,£| 


(2-31) 


Equation  (2-31)  is  a  curve- fit  to  the  law-of- the -wall  itself, 
Eqn. (2-15),  and  Ye  is  needed  to  evaluate  the  pressure  gradient 
parameter  Q  above.  Note  that,  in  the  limit  of  very  large 
r*  ,  Q  approaches  (a6+),  which  is  the  incompressible  two- 
dimensional  parameter  originally  used  by  White  (1969)* 
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The  parameter  "A"  in  these  formulas  accounts  for  the 
combined  effect  of  well  temperature  and  Mach  number,  as 
was  done  in  White  and  Christoph  (1972).  It  is  the  same 
parameter  defined  in  a  flat  plate  analysis  by  van  Eriest 
(1956): 

A  -  (Ta)/Te  -  1)  /fsin'1(a/c)  +  sin'-^b/c)!  ,  (2-38) 
where  a  =  <Taw«w)/re  -  2  ;  b  =  (T^-Tj/T,  . 

c2  .  r(TaH+Tw)/Te]£  -  4Tw/Te  . 

Values  of  A  for  various  Mach  numbers  and  wall  temperature 
ratios  are  shown  in  Figure  2-3*  Note  that  the  incompressible 
adiabatic  limit  is  A  =  1. 

2.4  Comparison  with  Experiment 

Equation  (2-25)  may  be  used  to  compute  the  turbulent 
skin  friction  distribution  C^(x)  =  2/\2(x)  along  a  body  of 
revolution  in  arbitrary  subsonic  or  supersonic  flow,  provid¬ 
ing  only  that  variations  in  wall  temperature  are  not  too 
great.  A  starting  value  X  (x=x  )  is  needed;  if  the 
computation  is  started  at  the  transition  point,  the  value 
=  20  is  recommended.  The  solution  of  (2-25)  may  be 
continued  downstream  until  the  separation  point  (if  any)  is 
predicted  to  occur  when  the  coefficient  (A.-3CXA2)  vanishes. 

We  now  consider  three  applications:  a)  flow  along  a 
slender  cylinder;  b)  supersonic  attached  flow  past  a  cone; 
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Figure  2.3.  THE  COMPRESSIBILITY  AND  WALL 

"TEMPERATURE  PARAMETER  A,  FROM  EQN.  (2-32)  . 
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and  c)  supersonic  -'lew  past  a  slender  body  of  revolution 


2.^.1  Compressible  Flow  Along  a  Slender  Cylinder 

The  application  of  the  present  theory  is  to  consider 

the  effect  of  compressibility  and  transverse  curvature  with 

no  pressure  gradient-  A  evlindsr  has  r  =  constant  and,  to 

o 

excellent  approximation,  Up  =  UQ  =  constant.  Equation  (2-25) 
reduces  to: 


A1  dx* 


*L 


(2-33) 


which  may  be  separated  and  integrated  for  a  direct,  relation 
between  skin  friction  and  local  Reynolds  number: 


(»Vue)(T 


>  + 
j*  A1(>,0,e,v,ro)  6\  (2-3^) 


With  zero  pressure  gradient,  the  curve-fit  approximation 
to  A1  is: 

A,  =  8.0  expj'j  (0.48  1  -  O.JO  Z  +  0.0018  Z‘)l  ,  (2-35) 

*• A  J 

where  Z  =  exp(0.4  VA)/r* 

Further,  rQ  is  directly  related  to  >  and  the  radius  Reynolds 

number  R  =  U  r  /v  ,  as  follows: 
a  e  o  e 

ro  * 

Since  R  is  constant  for  the  cylinder,  Eqrs.  (2-35,36) 
be  substituted  into  Eon.  (2-3*0  arid  the  integration  err; led 
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out  for  the  desired  relation  ^X(M»  or  better,  its  inverse 
C-.  =  C_(R  ,  R  ).  As  an  example,  the  computed  skin  friction 
on  an  adiabatic  cylinder  at  a  Mach  number  =  5  is  shown 
in  Figure  2-4.  It  is  seen  that  there  are  substantial 
transverse  curvature  effects  which  raise  the  skin  friction 
over  its  flat  plate  value.  As  R  becomes  very  large,  which 
is  equivalent  to  the  cylinder  being  very  short,  the  flat 
plate  skin  friction  is  approached  from  above. 

Since  the  integral  in  Eqn.  (2-34)  cannot  be  found  in 
closed  foi.i.  it  is  appropriate  to  find  a  curve-fit  formula 
for  the  results,  in  the  spirit  of  the  correlations  (2-1) 
proposed  for  incompressible  flow.  It  was  found  that  the 
flat  plate  formula  of  White  ruid  Christoph  (1972)  could  be 
gene  raided  to  Include  the  effect  of  x/rQ  =  Rx/Rft  as 
follows : 

1 

cf  =  0.455/[A2an2(b  R^  (V*w>  <Vt//A)  ",  (2-37) 

where  b  =  0.06/fl  +  0.025(x/r  )^’~|  . 

L-  °  j 

This  formula  is  valid  with  accuracy  to  +10y  over  the  entire 
turbulent.  Reynolds  number  range  (lCp-10-*)  and  for  Mach 
numbers  from  zero  to  ten.  The  slenderness  ratio  (x/rQ)  may 

r 

be  as  high  as  10°;  note  that  the  formula  predicts  a  roughly 
556  increase  in  skin  friction  for  (x/rQ)  «s  low  as  ten.  For 
very  small  (x/ro),  the  formula  reduces  to  the  compressible 
flat  plate  formula.  Equation  (24)  of  White  and  Christoph 
(1972) • 
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Equations  (2-  ?--)  or  (f-3?)  nay  be  compared  with 
experimental  data  of  Richmond  v!957)  for  supersonic 
turbulent  flow  past  long  cylinders.  Richmond  generated  a 
supersonic  flow  at  M  =  5-6  past  three  cylinders:  rQ  =  0.012, 
0.032,  and  0.125  inches,  respectively.  For  the  largest 
cylinder  (0.125 skin  friction  was  measured  with  a  floating 
element  balance  for  one  condition  (R  =  20,400)  and  was 
estimated  from  the  velocity  profiles  for  the  other  two 
cylinders.  The  authors  attempted  an  alternate  skin 
friction  determination,  using  a  sort  of  axisymnetric 
"Clauser-plot"  which  may  be  inferred  from  our  postulated 
law-of-the-wall,  Eqn.  (2-15) «  Since  the  cylinders  are 
assumed  adiabatic  and  the  pressure  gradient  is  zero  (a=(3=0), 
Eqn.  (2-15)  may  be  integrated  in  closed  form  to  give  the 
following  expression: 

u+(Y+,0,0,y,rp  -  £  sir"  *  m(Y+/Y^)'1  ,  (2-38) 

where  Y+  =  0.1108  ~o  ma^ch  with  the  incompressible  log  law 
c 

in  the  limit.  By  introducing  'he  definition  of  v  from 
Eqn.  (2-13)  ana  rearranging,  we  obtain  an  equivalent 
logarithmic  law: 


u/ur  =  i  ir(Y  /O.llt'8)  ,  where  u  =  sin  ^ (u/U^.) 


and  where  U  =  (2c  T  /r)^^ 
r  v  p  w  ' 


(2-39) 
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Figure  2-4.  SKIN  FRICTION  ON  AN  1  DIAS  ATI  C  CYLINDER  AT 


A  FREESTREAM  MACH  NUMBER  OF  5  FOR  VARIOUS  VALUES 
OF  R  =  U  r  /V  ,  COMPUTED  FROM  EQUATION  (2-34). 


k 


This  relation  shows  that  a  plot  of  the  "reduced"  velocity 

^  -A 

u  versus  Rao's  variable  Y‘  should  be  semi-logarithmic  and 
follo-f  the  incompressible  law.  Further,  if  the  velocity 
profile  u(y)  is  known,  one  may  use  Eqn.  (2-39)  to  infer  the 
friction  velocity  u*  and  hence  C^.  itself.  That  is,  if  we 
eliminate  u*  in  favor  of  Eqn.  (2-39)  niay  be  rewritten 
as  follows: 

C  =  /C^ 72  i  tn(r  STJZ)  +  So’  ,  (2-40) 

where  r  =  (3/Uj  (Tg/Tw)1//r  and  v  «  Rft  tn(r/rQ) 


Equation  (2-40)  may  be  plotted  with  as  a  par  neter,  as 
shown  in  Figure  2-5*  stfhen  data  for  C  versus  v  is  placed  on 
this  chart,  it  should  possess  a  logarithmic  overlap  layer 
from  which  one  can  infer  Cr  by  interpolating  between  the 
plotted  parametric  lines.  This  has  been  dene  in  Figure  2-5 
for  four  supersonic  cylindrical  profiles  given  by  Richmond 
(1957).  Note  that  all  profiles  demonstrate  the  proper 
slope  and  behavior  and  thus  substantiate  the  present  law-of- 
the-wall  epproach.  In  all  four  cases,  the  inferred  Cf  is 
lower  than  that  reported  by  Richmond  (1957).  Table  2.1 
summarizes  the  skin  friction  data  and  compares  with  two 
theories:  the  more  "exact"  integral  computation  of  Eqn.  (2-34) 
and  the  curve-fit  formula  of  Eqn.  (2-37).  Also  included  for 
comparison  is  a  flit  plate  computation  for  in  the  same 
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CLAUSER-PLOT  (EQ.  2-40). 


range-  Note  tnat  transverse  curvature  has  increased  C^. 

by  at  least  40$  to  as  much  as  100$.  The  present  theory 
seems  in  good  agreement,  falling  between  the  two  skin 
friction  estimates  (  Richmond  and  the  :'Clauser-plot'’ ) .  To 
the  author’s  knowledge,  no  other  theory,  whether  integral 
or  finite  difference  in  nature,  has  been  applied  to  these 
important  thick  axisymnetric  supersonic  boundary  layer 
experiments.  -Richmond  (1957)  also  reported  low  speed 
friction  and  velocity  profile  measurements  along  slender 
cylinders;  these  were  analyzed  in  the  incompressible  theory 
of  White  (1972)- 

2.4.2  Supersonic  Attached  Plow  Past  a  Cone 
A  classic  problem  in  the  literature  is  that  of  super¬ 
sonic  flow  past  a  cone  at  zero  incidence.  If  the  resulting 
shock  wave  is  attached  to  the  cone  vertex,  the  freestream 
flow  along  the  cone  surface  is  approximately  constant 
velocity  (a  =  0).  This  is  somewhat  analogous  to  the 
cylinder  case  of  the  previous  example,  except  that  here  the 

surface  radius  r  varies.  If  x  is  along  the  cone  generators, 
o 

we  have 


r^  =  x  sin  ®  ,  (2-41) 

where  ©  is  the  cone  half-angle.  It  is  of  interest  to  try 
and  relate  the  cone  skin  friction  to  an  equivalent  flat 
plate  c~  cylinder  at  the  same  Reynolds  number  R^.  Accordi-.r 
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to  the  classic  turbulent  analysis  by  van  Driest.  (1952), 
at  position  R^  on  the  cone  is  equal  to  on  a  flat 
plate  at  position  (R^/2).  Thus,  for  equal  R^,  the  cone 
skin  friction  is  10-15%  higher  than  the  plate  C?.  This 
analysis  of  course  does  not  account  for  boundary  layer 
thickness  effects,  which  the  present  theory  includes  as 
the  (r*)  effect  cn  the  law-of- the -wall.  With  U  =  constant, 
Eqn.  (2-25)  reduces  for  cone  flow  to: 


d\  -•  dro 

ft  +  a  ~  _ L 

A1  6x*  3  rQ  dx* 


(2-12) 


and,  since  is  constant,  (l/r0)  (drQ/dx*)  =  1/x*  from 
Eqn.  (2-11).  If  transverse  curvature  is  important,  i.e. 
if  the  boundary  layer  is  thick  compared  to  the  cone  radius, 
then  and  A^  depend  upon  r^(x,x*)  and  the  variables  are 
not  separable  in  Eqn.  (2-12).  However,  if  we  assume  a 
thin  boundary  layer  (large  r*),  then  our  curve-fit 
approximations  give: 

A^  -  ^  A,  =  a  e^  (a  =  £.0,  b  =C.l£/A) 

(2-43) 


for  this  case,  Eqn.  (2-12)  has  a  closed  form  solution: 


(P-lu) 


which  may  be  compared  to  the  analogous  solution  for  flat 
plate  flow,  Eqn.  (2-37)  or  (1-36): 
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X 


1 

b 


(2-46) 


in. 

D 


RjXVaJ 


By  comparing  (2-44)  and  (2-45),  one  may  deduce  a  first- 
order  turbulent  c;>ne  rule": 

C„  (at  Rj  =  cr  (at  IL/l  1+0*96/ A*)  (2-46) 

Icone  x  1  plate  *  ~ 

This  may  be  compared  with  the  classic  cone  rule  developed 
by  van  Driest  (1952): 

cf  (Rx)  =  c  (  V2)  C2-47) 

cone  plate 

From  Figure  2-3,  A  =  1.0  at  low  Mach  numbers,  so  that  the 
correction  factor  for  Reynolds  number  equals  1. 96  approxi¬ 
mately,  or  very  close  to  van  Driest* s  factor  of  2.0.  A 
later  analysis  by  Tetervin  (1969)  gives  a  factor  of 
(2.0  +  1/N),  where  N  is  the  exponent  in  an  assumed  (l/Nt  ) 
power-law  relation  between  C ^  and  momentum  thickness 
Reynolds  number.  Since  N  is  of  the  order  of  5  to  9» 
increasing  with  Reynolds  number  and  Mach  number,  Tetervin* s 
factor  is  also  abcut  two.  Figure  2-6. a  shows  values  of 
the  ratio  of  cone  to  plate  friction  coefficient  as 
calculated  from  Eqn.  (2-46).  The  predicted  increase  in 
cone  friction  is  of  the  order  of  5  to  15&  being  higher 
at  lew  Reynolds  numbers. 

In  fact,  however,  neither  Eqn.  (2-46)  nor  the  van 
Driest  or  Tetervin  analysis  is  accurate  for  small  cone 
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a)  THIN  BOUNDARY  LAYER  ASSUMPTION:  EQN, (2-46). 


Figure  2.6.  ILLUSTRATION  OF  THE  TURBULENT  CONE  RULE. 
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angles,  because  of  the  transverse  curvature  effect 
(large  6/rQ).  That  is,  a  numerical  or  graphical  solution 
of  Eqn. (2-42) ,  using  the  actual  formulas  fcr  and  from 
Eqns. (2-29, 31),  shows  a  further  increase  in  cone  skin  friction 
at  the  smaller  cone  angles.  These  numerical  solutions  are 
illustrated  in  Figure  2.6.b  for  a  freestream  Mach  number  of 
2.0  along  the  cone  surface.  Thus  transverse  curvature  adds 
an  additional  10  to  30#  to  the  skin  friction,  and  only  the 
thicker  cones  (20°  to  45°  half  angle)  approximate  the  first 
order  cone  rule  above  (2-46,47).  By  analogy  with  the  cylinder 
analysis  (Eqn. 2- 37),  the  relevant  transverse  curvature 
parameter  is  (x/rQ)  =  csc(^)  for  a  cone.  Hence  the  numerical 
solutions  in  Figure  2.6.b  may  be  correlate  into  the  following 
formula  for  constant  pressure  (attached)  flow  along  a  cone: 


:f(c.one)  = 


0.455 

>  2  n  !  in  At 


A"  lb  (ue/uw)/A 


1 


(2-48) 


where  b  =  0.06(1  +  0.96/A)*1(l  +  0.03  esc  y)”1, 


This  formula  i~  accurate  to  about  +3^  over  the  complete 
supersonic  turbulent,  cone  flow  range,  including  constant 
nonsdiabatic  wail  temperature,  which  is  accounted  for  by 
the  van  Driest  parameter  A  in  the  formula  (Eqn. 2-32). 
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TABLE  2- 


2.4.3  Supersonic  Flow  with  Pressure  Gradient: 

For  flow  past  a  general  body  of  revolution,  the  complete 

differential  equation  for  skin  friction,  Eqn.(2-25),  should 

be  used,  since  a  t  0  and  r  =  r  (x).  In  one  sense,  these 

o  o 

effects  are  mutually  exclusive,  because  if  pressure  gradients 
are  important,  transverse  curvature  is  probably  not  important, 
and  vice  versa.  That  is,  only  a  slender  body  can  have  large 
(6/ro),  but  a  slender  body  cannot  have  large  pressure  gradients. 
However,  the  shape  effect  term  (A^/ro) (drQ/dx*)  is  nearly 
always  important. 

We  give  here  a  typical  solution  of  the  general  theory 
by  comparing  Eqn.(2-25)  with  the  skin  friction  measurements 
of  Allen  (1970)  for  supersonic  flow  past  a  so-called  Haack- 
Adams  body  of  revolution,  shown  to  scale  in  Figure  2.7.  The 
shape  of  the  body  is  given  by  the  following  formula: 


r  /r 
o  max 


O.TCPtl-t2)1*5*  0.16934  cos“At  +  t(l-t2)1/2 


where  t  =  (1  -  2Z/L) . 


(2-49) 


Note  that  Z  Is  the  axial,  not  the  surface,  coordinate  (see 

Figure  2.1).  Allen  (1970)  tested  a  wjnd  tunnel  model  with 

L  =  36  inches  and  r  -  l.S  inches.  Skin  friction  at  seven 

max 

stations  was  measured  with  a  Preston  tube  and  also  estimated 
by  1)  the  3aronti-Libby  law;  and  2)  the  Fent^r-Stalmach  law. 
The  upstream  Mach,  number  was.  varied  from  2.5  to  4.5,  and 


the  example  shown  (M^  =  2.?6)  is  typical  of  both  the  data 
and  the  theory.  With  L  =  36  inches,  the  Reynolds  number  RT 
in  Eqn. (2-25)  was  equal  to  3*0x10^.  The  freestream  Mach 
number  distribution  Mg(x)  is  shown,  and  the  velocity 
distribution  is  of  a  similar  shape  and  was  curve-fit  by 
the  expression 

V(x*j  =  1  +  0.156(t  -  t2)  ,  t  =  x*  -  0.139*  (2-50) 

Figure  2.1  shows  three  different  analytical  estimates: 

1)  the  exact  theory,  Eqn. (2-25);  2)  a  flat  plate  theory 
based  on  local  R^,  Eqn. (I-36) ;  and  3)  a  long  cylinder 
theory,  Eqn. (2-37),  based  on  local  R  and  r  (x).  It  is 
seen  that  the  cylinder  and  flat  plate  formulas  give  an 
estimate  of  the  average  skin  friction  on  the  body.  The 
difference  between  these  two  is  about  5#>  which  represents 
the  transverse  curvature  effect.  The  exact  theory  is  in 
excellent  agreement  with  the  data  over  the  entire  body. 

The  higher  values  of  C^.  (compared  to  a  flat  plate)  at  the 
front  of  the  body  are  due  to  roughly  equal  contributions 
from  1)  the  term  (A^/rQ) (drQ/dx*) ;  and  2)  the  favorable 
pressure  gradient.  Similarly,  the  lower  values  at  the  rear 
are  due  to  1)  a  change  in  sign  of  (drQ/dx*);  and  2)  the 
adverse  pressure  gradient.  The  same  excellent  agreement 
with  experiment  was  found  for  the  other  three  conditions 
(M  =  2.5,  3*95,  4.5)  investigated  by  Allen  (1970). 
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DATA,  ALLEN  (1970) 


\>  0 


o  o 


o 


o 


0_J _ | _ I _ I _ L 

0.2  0.4  0.6  0.8  1.0 

x*=  x/L 


DATA,  ALLEN  (1970): 


M  =2.96 


0.2  0.4  0.6  0.8  1.0 

x*=  x/L 


Figure  2.7.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 

SUPERSONIC  FLOW  PAST  A  BODY  OF  REVOLUTION  AFTER  ALLEN  (1970) 


2.5  Summary 

This  chapter  has  presented  a  complete  theory  for  the 
computation  of  turbulent  skin  friction  in  compressible 
axisymmetric  flow  past  arbitrary  bodies  of  revolution.  The 
approach  leads  to  a  single  first  order  ordinary  differential 
equation,  Eqn.(2-25)>  for  Cf.(;  )  and  accounts  for  supersonic 

i 

flow,  nonadiabatic  wall  temperatuie,  and  extremely  thick 
boundary  layers,  where  a  transverse  curvature  correction  is 
necessary.  Applications  to  cylinder  flow,  cone  flow,  and 
flow  past  a  pointed  body  of  revolution  all  show  good 
agreement  with  experiment.  It  is  therefore  thought  that 
the  present  theory  has  no  significant  deficiencies.  Since 
this  theory  is  also  apparently  the  simplest  comparable 
analysis  to  be  found  in  the  literature,  we  therefore 
recommend  it  to  engineers  for  general  usage. 
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Chapter  Three 


The  Prediction  of  Three  Dimensional  Skin  Friction 

3-1  Introduction 

The  most  common  type  of  boundary  layer  encountered  in 
practice,  and  the  most  intractable  from  the  standpoint  of 
analysis,  is  the  three-dimensional  turbulent  boundary  layer. 
Reviews  of  the  state  of  the  art  with  respect  to  this  impor¬ 
tant  class  of  flows  have  been  given  by  Cocke  and  Hall  (1962), 
Joubert,  Perry  and  Brown  (1967)  and  more  recently  by  Nash 
and  Patel  (1972).  Existing  theories  are  few  and  seem  unduly 
complicated  while  available  experimental  data  are  usually 
unreliable.  Historically,  two  approaches  have  been  used  to 
analyze  three  dimensional  turbulent  boundary  layers;  the 
integral  approach  and  the  differential  approach. 

Typical  of  the  integral  approach  is  the  work  of  Smith 
(1966),  and  Cumpsty  and  Head  (1967)-  These  authors  used  two 
momentum  integral  equations  plus  an  auxiliary  equation  to 
account  for  variations  in  the  st Teamwise  shape  factor.  They 
extended  Head’s  (1958)  two-dimensional  entrainment  relation 
to  three  dimensions  but  retained  his  two-dimensional  empirical 
functions. 

Differential  methods  have  been  given  by  Bradshaw  (1971) 
and  Nash  (1969);  both  are  extensions  of  the  two-dimensional 
method  of  Bradshaw  et  al  (19^7)  which  utilizes  the  turbulent 
energy  equation.  The  major  difficulty  with  differential 
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methods  is  in  the  specification  of  boundary  conditions.  For 
a  given  domain,  boundary  conditions  must  be  specified  along 
all  surfaces  where  fluid  is  entering  from  the  outside. 

In  this  chapter  a  new  integral  method  is  developed 
which  greatly  reduces  the  computational  difficulties  inherent 
in  the  calculation  of  three  dimensional  skin  friction.  The 
analysis  results  in  two  coupled,  non  linear  partial  differ¬ 
ential  equations  with  the  skin  friction  coefficient  and  the 
tangent  of  the  angle  between  the  total  surface  shear  stress 
vector  and  the  shear  stress  vector  in  the  freestream  direction 
as  the  only  unknowns.  This  is  accomplished  by  assuming 
suitable  "law-of-the-wall"  velocity  correlations  for  both 
the  freestream  and  cross  flow  directions  end  then  inte¬ 
grating  the  governing  equations  with  respect  to  the  law-of- 
the-wall  coordinate  y+.  The  resulting  partial  differential 
equations  may  be  easily  solved  by  finite  difference  tech¬ 
niques  since  they  require  initial  data  along  only  a  single 
curve. 

Thus  the  new  method  Is  a  compromise  between  the  "classi¬ 
cal”  integral  and  differential  techniques.  It  is  much  more 
straight-forward  and  contains  considerably  less  empirical 
content  than  the  former  and  is  computationally  much 
simpler  than  the  latter. 

At  the  present  time  this  method  is  limited  to  incompress¬ 
ible  flows;  however,  extending  it  to  account  for  compress¬ 
ibility,  heat  transfer,  and  variable  fluid  properties  should 
be  possible  in  view  of  chapter  1  of  this  report  and  the  past 
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efforts  of  White  and  Christoph  (1970). 

In  order  to  discuss  three-dimensional  boundary  layer 
flows  it  is  most  appropriate  to  work  In  a  curvilinear 
coordinate  system  composed,  in  part,  of  the  projection  of 
the  free  streamlines  onto  the  boundary  surface.  The 
coordinate  along  chese  free  streamline  curves  is  designat¬ 
ed  as  s,  while  the  coordinate  along  the  orthogonal 
trajectories  of  the  free  streamlines,  in  the  surface,  is 
designated  as  n,  and  the  coordinate  normal  to  the  surface 
is  designated  as  y.  Corresponding  to  these  three  coordinate 
directions  are  the  metric  coefficients  h-^,  hg,  and  h^,  and 
the  mean  velocity  components  u,  w,  and  v  respectively.  For 
simplicity,  it  shall  be  assumed  that  the  curvature  of  the 
bounding  surface  does  not  change  abruptly  and  that  the 
boundary  layer  thickness  is  small  compared  with  the  principal 
radii  of  curvature  of  the  bounding  surface,  so  that  h^  may 
be  taken  equal  to  unity  and  y  becomes  a  simple  distance 
normal  to  the  surface.  As  a  consequence,  h^  and  hg  are 
known  functions  of  s  and  n  provided  that  the  inviscid  flow 
over  the  surface  is  known.  This  coordinate  system  is  shown 
in  Figure  J . I . 


1  2l- 
hjlxj  >n 


the  s  component  of  the  momentum  equation 


(3-1) 
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(3-2) 


u  ~ji  +  w  ~ii  +  v  ~  +  uy  _  wg  ahg 

*s  EJ  >n  *y  h^'  ~  HpC]’  *s 

=  .  _i_*£  +  1  Us 

ph1  *s  p  *y 

and  the  n  component  of  the  momentum  equation 

U  vw  +  W_  >Jrf  +  aw  uw  an2  _  u  1 

E|~  *s  h?  *n  h-jh^  as  Kjlv  >n 

=  .  4_  »E  +  it2i  •  (3-3) 

ph2  *n  p  *y 

In  these  equations  Ts  and  Tn  represent  total  shearing 
stresses  and  induce  the  effects  of  viscosity  and  Reynolds 
stress.  Since  only  incompressible  flow  is  to  be  considered, 
the  energy  equation  is  not  needed. 

The  distinguishing  feature  of  a  three  dimensional 
boundary  layer  is  crossflow  (the  w  component  of  velocity) 
and  the  resulting  vector  character  of  the  surface  shear 
stress.  Crossflow  occurs  when  the  frees tream  streamlines 
are  curved,  which  gives  rise  to  an  unbalanced  centrifugal 
pressure  gradient  in  the  boundary  layer.  Such  crossflow 
imposes  a  shear  stress  on  the  surface,  normal  to  the  free- 
stream  direction  (T  ),  which  skews  the  total  surface  stress 
vector.  If  the  direction  of  curvature  of  the  frees tream 
streamline  is  reversed,  skewing  in  two  lateral  directions 
can  occur  in  a  single  velocity  profile.  Typical  velocity 
profiles  showing  skewing  in  one  and  two  directions  are  shown 
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Figure  3.2,.  THREE-DIMENSIONAL  BOUNDARY  LAYER 


WITH  BILATERAL  SKEWING. 


in  Figures  3-1  and  3*2 • 

Since  the  present  method  depends  upon  an  apriori 
specifica"ion  of  velocity  correlations  in  both  the  free- 
stream  direction  and  the  crossflow  direction  it  would  seem 
appropriate  to  discuss  these  in  some  detail. 


?.?  Specification  of  Velocity  Correlations 

In  general,  integral  techniques  are  insensitive  to 
derails  of  the  assumed  velocity  correlations.  However, 
this  does  not  mean  that  one  can  be  completely  cavalier  in 
their  specification-  In  order  to  insure  that  the  chosen 
correlations  incorporate  the  correct  physical  parameters, 
consider  that  near  the  wall  the  force  balances  given  by 
equations  (3-2)  and  (3-3)  reduce  to 


and 


(3-.U) 


(3-5) 


Assuming  that  the  shear  stress  in  the  freestream 
direction  in  a  three-dimensional  boundary  layer  may  be 
represented  by  the  same  eddy  viscosity  relation  as  for  a 
two  dimensional  boundary  layer  gives 

Ts  =  p,c2y2l|y5|7  (3-6) 

with  von  Kerman's  constant  k  equal  to  0.4.  Nondimens ional- 
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izing  equations  (3-4)  and  (3-6)  with  respect  to  wall 

quantities  and  solving  for  the  velocity  gradient  gives 

+  +  1/2 
au.  _  (1  +  ay  ) 

W  ~  - ~- 

K  y 


The  pressure  gradient  parameter  is  defined  by 


a 

s 


If 


(3-E) 


1  /O 

where  u*  =  (Tsw/p)  '  is  the  shear  velocity  in  the  free— 

stream  tiiiection.  Equations  (3-7)  and  (3-8)  are  exactly 

the  same  results  as  have  been  previously  obtained  by  White 

and  Christoph  (1972)  for  the  two-dimensional  case;  except 

that  here  a  includes  the  curvilinear  scale  factor  h, . 

s  1 

Therefore,  it  will  be  assumed  that  the  law-of-the-wall 
correlation  for  the  streamwise  velc.'ity  component,  in  a 
three-dimensional  boundary  layer  (Eqn.  ~-7)  does  not  differ 
from  the  two-dimensional  correlation  r,n*is  is  a  reasonable 
assumption  as  has  been  pointed  out  K,  .  ash  (1969)  and  others. 

It  would  bp  extremely  attractive  to  carry  through 
similar  arguments  for  the  c-ossflow  correlation.  However, 
this  immediately  leads  to  excessive  complication.  The 
reason  is  that  the  simple  eddy  viscosity  approximation  of 
equation  (3-6)  does  not  appear  adequate  to  describe  the 
crossflow  profile.  Instead  one  must  invoke  a  coupled  eddy 
viscosity  expression  such  as  suggested  by  Cousteix,  vjuemard 
and  Michel  (1971); 
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2,2.au.dw 


n 


!dy'^y 


(3-=?. 


where  F  is  a  sublayer  damping  factor  suggested  by  van  Driest 
(1956) 


=  1-6X13  [_25qT  (TP)x/2j 


(3-10'i 


end 


1  -  0.085  6  Tanh  (575^  *^p)  (3-11; 


Nondimensionaiizing  and  combining  (3-5)  and  (3-9) 
leads  to 


>  +  ,  4-  * 

=  (»  *  an  y  > 

liSll  iV 
>y  ! 


(3-^1 


with 


a  =  — 

n  t 


1 


IA 

T1W  2 


(3-1:5 


The  symbol  0  stands  for  the  tangent  of  the  angle  between 
the  total  surface  shear  vector  and  the  shear  stress  vector 
in  the  freestream  direction. 

Expression  (3-1?)  is  obvi  'usly  much  too  complicated  tr-. 
be  tractable  as  a  crossflow  velocity  correlation.  Fortuna* -.V 
an  alternative  exists  in  the  simple  hodogrsph  models,  which 
often  give  a  better  representation  of  the  crossflow  velocity 
profile  than  does  (3-12). 

The  earliest  hodograph  model  was  proposed  by  Prandtl 
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(19^*6)  who  speculated  that 


w  =  ueg(|)  (3-14) 

where  g  was  a  function  such  that  g(0)  =  1  and  g(l)  =  0.  He 
took  g  to  he 

«(f)  =  1-  $  (3-15) 

Mager  (1951)  suggested  a  more  accurate  form  in  which 

g(f)  =  (1-  (3-16) 

Recently,  several  more  sophisticated  hodograph  models  have 
been  proposed.  Johnston  (1957)  represented  the  crossflow 
velocity  by  a  triangular  hodograph  such  that 


u 

U 


<ir> 

e  p 


(3-17) 


(3- 


where  (u/D  )  corresponds  to  the  apex  of  the  triangle  and  is 

G  p 

given  by 


(1  + 


(3-18) 


and  U  is  the  freestream  velocity.  Johnston's  model  is 
sketched  in  Figure  3-3-  This  model  has  several  shortcomings. 
There  are  situations  in  which  the  outer  region  is  not 
adequately  represented  by  a  straight,  line  and  there  are 
situations  ir.  which  the  apex  of  the  triangle  is  not  well 
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Figure  3.3.  SKETCH  OF  THE  TRIANGULAR  HOLOGRAPH 


APPROXIMATION,  AFTER  JOHNSTON  (3957) 


defined.  Also,  he  experimental  determination  of  P  is 
very  difficult  due  to  the  fact  that  the  inner  region  of  the 
model  actually  corresponds  to  the  viscous  sub  layer. 
Eichelbrenner  and  Peube  (1966)  modeled  the  crossflow  velocity 
by  a  polynomial  of  the  form 


w  u  .  Ft  .  /U  x 

FT®  1  +  a(u~' 


y  f  U  /  ,U 

b'TJ")  +  cvu“) 
c  e 


A  ,  u  %  *  ,  ,  a  x 

d(y  )  +  e ( (j  ) 

e  e 


5-i 


(3-19) 


where  a,b,c,d,e  are  evaluated  by  using  boundary  conditions 
at  the  wall  and  at  the  outer  edge  of  the  boundary  layer. 
These  constants  are  also  extremely  difficult  to  verify 
experimentally.  This  model  does,  however,  have  the 
advantage  of  being  able  to  predict  the  "S-shaped”  crossflows 
which  discussed  in  Section  3-1  and  observed  experi¬ 

mentally  by  KLinksIek  and  Pierce  (19^0). 

With  these  considerations  in  mind  it  was  decided  by 
the  authors  to  employ  the  hodograph  model  as  proposed  by 
Mager  (1951)  in  the  form 


-I 


(3-20) 


This  choice  seemed  a  good  compromise  between  complexity  and 
relative  accuracy.  In  fact  this  correlation  compares  well 
with  some  experimental  data  as  shown  in  Figure  3-5-  Until 
such  time  as  more  reliable  three-dimensional  data  become 
available,  and  in  particular  profile  data,  the  use  of  a  more 
sophisticated  hodograph  model  does  not  seem  warranted. 
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-•3  The  Method  of  Solution 

In  order  to  solve  equations  (3-1),  (3-2),  and  (3-3) 
they  are  first  nondimensionalized  using  the  shear  velocity 
in  the  freestream  direction.  The  result  is 


u*u+  a 


•}* 

*U  3  /  w.  +v  ,  U*w  3  /  -  +  %  ,  U*  t 

Ts  (u  *  1 ZT  dn  *u  u  )  v  v“  (u  u  ) 

t  -£  uV  lh  u*V2  ah;  _  _  1  aP  .  U*  aTs 

u  TTTEJ  3  n  “  h  hu  rs  ph,  Ts  pv  -'V f 

1  1  x 


u*  d 


(?-?l) 


and 


u*u+  a  ,  -  +,  u*w+  a  ,  _  +,  u*  a 

7 -5?  <u  a  )  "h^"  1Sn  (a  ^  )  +  v  “  ^  <U/?  ) 
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u*?  u+w+  _^2  _  u*V2  ahl 


1  dp  u*  ®Tn 
"Sn  +  nv  ‘5Tr+ 
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where  the  velocity  v  is  determined  from  the  continuity 
equation: 


v  f  r  1  3  1  a  ._  _  +."»  + 

v  =  '  u *  5  57  u  )  +  Kfe  TSn  (hlu^  > 


(3-23) 


Expressions  for  the  partial  derivatives  with  respect  to  both 
s  and  n  are  needed.  Note  that 

u+  =  fcn(y+,  ag)  (3-24) 

and 

w+  =  fcn(y+,  as,  6,  6+)  .  (3-25) 


Therefore, 
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XI  ^  + 
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as  dy 


da  - 
,  s  dw 

+  3s~35“ 
s 


d0  dvr+  ?6+dw1" 
+  ?¥  dS^ 


(3-26) 

(3-27) 


plus  similar  expressions  for  the  n  derivatives.  For  the 
purpose  of  skin  friction  calculations,  it  is  convenient  to 
make  the  following  non-dimensionalizations: 


(2/cf6)1/2j  h*  =  ly-L,  h-  .  tv/L,  rl  =  -2 


U  L 


(2-28) 


V  =  Ue(s,n)/Uc 


where  L  is  a  reference  length,  and  UQ  is  a  reference  speed. 
Then,  from  the  definition  of  a 

s 


+ 


(i/v) 


(3-29) 


and 


>a  3a  l-i  3  >2 

*n~  ~  ~ an-  +  H£~  (1/v> 


(3-30) 


The  partial  derivatives  of  6+  with  respect  to  s  and  n  are 
more  difficult  to  compute.  Equation  (3-7)  Integrates 
analytically  to 


=  ~  r^«p-i)  +  <Ur  (sV» 

SO 


(3-31? 
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where 


p  =  (i  +  asy+)1/2 


At  the  edge  of  the  boundary  layer,  equation  (3-31)  becomes 

u+  .  x  =  1  r2(p  -  l)  +  ln(p6-^  (_V))1  (3-32) 

e  1  *  -  6  VT  - 


where 


=  (1  +  “  6+) 


1/2 


Differentiating  equation  (3-32)  with  respect  to  s  and  n 
results  in: 
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(3-33) 


and  a  similar  expression  for  the  n  derivative.  Combining 
equations  (5-26),  (3-27),  (3-29),  (3-30),  and  (3-33)  gives: 
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After  substituting  the  partial  derivatives,  equations  (3-3*0 
through  (3-37),  into  equations  (3-21)  and  (3-22)  and  inte- 
nrntir.-  the  resulting  equations  from  y+  =■  0  (7  =  t  )  to 
y+  =  f>+  (t  -  0).  one  obtains 


78 


Wt*.***1 


1  SV  2„  + 


1  >V 


%  ?T  (Gl-3“s^>  +  v%  <X1  8  Gl)  -  7K*  5K  g3 

A 

,  3X,  X,  G„  2 

+  h*  •Sir  (G3-3a -'iii+'lG5^  ■  ^i-2  W 

xi4q^  *2  ,i\  *  yi  3£  -  _  h  >hl  G 

RTh*  >s>n  '  7'  hlf  >n  G6  Hfhf  ?n  G7 

L  c  X  d 


X1  ah2 


hfe  >r  Ge  =  V 


2  1 


(3-38) 


1  axl 


H*  ?T  ^G3"3asG9"G10)  4  W*  Is  (X1£?6+'G11)  -W*  !?  °3 


ki  ay 


H*  3~  ^Gll~^asG12”xlG13^  -  ^nf  It2  <4> 


-  4r 

Ai  Gi2  a 


(3-39) 


1  A  O 


j.  j.c  v  /j-s  1  a@  ,, 

RlS|“  "Sn^s  ^V1  "  SJ  a?  Gl4  "  h|  an  G15 


'1  ah! 


x1  ah* 


+  K*nf  aH“  gi6  "  S|EJ  —  °17  "  eRLV 


The  G  coefficients  in  equations  (3-38)  and  (3-39)  are 
defined  .in  terns  of  quadratures  over  the  velocity  correlations 
and  various  gradients  of  th*'  velocity  correlations.  These 
are  tabulated  and  evaluated  in  appendix  A. 

Equations  (3-38)  and  (3-39)  are  two,  coupled,  non¬ 
linear,  partial  differential  equations  which  have  for 
dependent  variables  only  the  skin  friction  coefficient  and 
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the  tangent  of  the  angle  between  the  total  surface  shear 

stress  vector  and  the  shear  stress  in  the  freestream 

direction  -  no  integral  thicknesses  need  be  calculated. 

The  inviscid  flow  solution  is  assumed  known,  a  may  be 

s 

computed  from  its  definition,  and  6  may  be  computed  from 

equation  (3-32).  As  should  be  the  case,  equation  (3-38) 

reduces  to  the  incompressible,  two-dimensional  analysis 

of  White  (1969)  when  =  0,  w+  =  0,  and  . 

2  1 

3.4  Comparison  with  Experiment 

3 . 4 . 3  The  Rotating  Disk  Problem 

A  convenient  starting  point  to  test  the  validity  of 
any  three-dimensional  theory  is  the  problem  of  a  smooth 
plane  disk  rotating  in  an  otherwise  stationary  fluid. 

While  axisymmetric,  this  problem  exhibits  the  essential 
feature  of  a  three-dimensional  flow,  namely  crossflow  and 
th«“  resulting  ve^or  character  of  the  surface  shear  stress. 
The  present  theory  has  been  applied  to  the  rotating  disk 
problem  in  a  paper  by  Lessraann  and  Christoph  (1972)  and 
is  discussed  below. 

The  boundary  layer  equations,  equations  (3~1)>  (3-2) 
and  (3-3),  may  be  reduced  to  those  governing  the  boundary 
layer  flow  on  a  rotating  disk  by  identifying  s  with  the 
coordinate  and  n  with  the  coordinate  r,  both  fixed  in 
space.  The  metric  coefficients  become  h^  =  r  and  h^  =  1 
and  the  equations  reduce  to 
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(3-40) 


and 
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All  derivatives  with  respect  to  0^  have  been  eliminated 
because  of  symmetry. 

The  boundary  conditions  on  the  equations  are  that  at 
the  rotating  plate  v=w=0  and  u=-rf.  Here  0  is  the  angular 
velocity  of  the  rotating  disk,  assumed  counterclockwise 
about  the  y  axis.  In  the  freestream  u,  v,  and  w  as  well  as 
the  shearing  stresses  are  zero.  As  a  consequence  of  this 


T?  = 


0 


(3-43) 


For  what  follows,  it  will  be  convenient  to  transform 
equations  (3-41)  and  (3-42)  so  as  to  obtain  more  conven¬ 
tional  boundary  conditions.  Defining  a  new  circumferential 
velocity  as 

u'  =  u  +  rO  (3-44) 

the  boundary  conditions  at  the  piste  become  u’  =  v=w=0  and 
in  the  freestream  v=w=0  while  u1  =  rf>.  Using  equation  (3-44) 
equations  (°-4l)  and  (3-42)  become 
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(3-45) 
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and 

*  I?  +  I?  +  2°“*  -  ***  =  ^  •  (3-46) 

Mager1 s  crossflow  velocity  correlation,  equation  (3-20), 

will  be  used  and  since  a  =0, 

s 

u+  =  2.51n(y+)  *  5o  (3-47) 

will  be  used  for  the  streamwise  velocity.  Comparison  of 
these  two  velocity  correlations  with  the  data  of  Cham  and 
Head  ( 1969)  is  shown  in  Figures  3-4  and  3-5*  The  velocity 
v  is  eliminated  from  equations  (3-45)  and  (3-46)  by  the 
continuity  equation,  and  the  resulting  equations  are  non- 
dimensionalized  with  respect  to  wall  variables.  For 
convenience  u*  is  related  to  the  circumferential  skin 
friction  coefficient  by 

X1  '  -  $  (3-«) 

and  r  is  non-dimensionalized  by  the  relation 

*  =  r(n/v)i,/2  .  (3-49) 

Note  that  v.  is  the  square  root  of  the  radial  Reynolds 

number.  In  integrating  the  two  momentum  equations  from 

y4"  =  0  to  y+  =  6^,  the  shear  velocJ  ty  in  the  radial 

1/2 

direction,  v*  =  '  •  will  appear.  The  shear  velocity 
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This  problem  reduces  to  the  consideration  of  coupled 
ordinary  differential  equations  as  a  result  of  the  fact  that 
the  rotating  disk  is  a  twc-cocrcir.ate  problem. 

The  results  of  the  present  analysis  consist  cf  pre¬ 
dictions,  using  equations  (3-51)  and  (3-5?)»  for  the 
circumferential  skin  friction  coefficient  c „  ,  and  the  ratio 
of  the  crossflow  shearing  stress  to  the  circunferentia’ 
shearing  stress,  S.  These  are  cospared  with  experimental 
results  ir.  Figures  3-6  and  3-7-  Figure  3-6  shows  that  the 
present  calculation  of  cr_  compares  very  favorably  with  the 
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experimental  cats  and  falls  approximately  halfway  between 
the  results  of  Cooper  (1971 }  and  von  Karman  (1946). 

Cooper  used  a  finite  difference  scheme  based  on  the  work 
of  Smith  and  "ebeci  (19^6 ,1971 )  *  Figure  3-7  is  reproduced 
rro:n  the  work  of  Cham  and  Head  (1969)  and  compares  several 
predictions  for  c  with  one  measured  value.  Unfortunately, 

*he  uncertainty  associated  with  this  single  data  point  is 
relatively  large,  k'ith  this  in  rair.d,  it  can  be  seen  that 
the  present  theory,  as  well  as  the  entrainment  method  of 
Cham  anu  Head  (1969),  and  the  prediction  of  Banks  and 
Gadd  (1962)  give  plausible  results.  Both  the  present 
calculation  and  that  of  Cham  and  Head  show  a  tendency  tc 
fall  with  increasing  radial  Reynolds  number,  while  the 
prediction  of  Banks  and  Gadd  (1962 )  is  for  a  constant  value 
of  ®.  Cooper  \lo~i )  din  not  concern  himself  with  predicting 
total  shear  and  hence  no  comparison  with  his  work  is  possib^^ 
on  this  polo.. 

Since  the  present  approach  and  the  analyst ^  of  von 
Kerman  {19^6 )  ere  integral  methods  it  is  informative  to 
consider  both  the  simi-srities  and  important  differences 
between  these  ~wo  works.  The  velocity  correlations  used 
by  von  Karman  are  (in  the  present  notation) 

1/7 

u  =  rn(v/6)  .  (3-59) 


and 


w  -  6u(l-v/5) 


(3-60) 


c9 


These  are  tc  be  compared  with  equations  (5~47)  ana  (3-20). 

It  is  well  known  that  the  1/7  power  law  gives  a  good 

representation  of  fully  turbulent  /elocity  profiles;  at 

least,  for  a  limited  range  of  Reynolds  numoer.  This  fact 

accounts  for  the  relatively  good  agreement  with  the 

circumferentiaj.  skin  friction  date.  At  the  same  time,  the 

crossflow  correlation  used  by  wn  Karman  does  not  give  e 

good  representation  of  the  radial  velocity  profile  as  is 

shown  in  Figure  (3-8)*  Also,  in  order  to  get  an  exact 

solution  of  his  equations,  von  Kaiman  assumed  p  to  be 

constant  and  the  boundary  layer  thickness  5  proportional  to 

r~  -  These  facts  account  for  the  relatively  poor  pre~ 

diction  of  9  =  0.162.  In  order  to  carry  through  the 

calculation  following  von  Karman,  one  must  express  the  wall 

shear  stresses  t.  and  t  as  functions  of  p  and  6.  This 
8w  rw 

is  in  contrast  to  the  piesent  analysis  in  which  the 
dependent  variables  are  9  and  the  circuafe  ’-ential  skin 
friction  coefficient.  No  shear  stress  correlations  are 
necessary.  One  can  carry  cut  von  Kerman's  analysis  retain¬ 
ing  the  r  dependence  of  £  and  6  and  also  consider  the  effect 
cf  replacing  (3-60)  with  Mager's  correlation.  As  is  shown 
in  Figure  3-6,  both  of  these  lead  to  essentially  the  same 
skin  friction  curve,  which  is  somewhat  low.  The  effects  of 
these  changes  on  the  prediction  of  9  is  more  drastic  as 
shown  in  Figure  3-7.  When  one  keeps  the  original  velocity 
correlation  but  allows  fo?'  r  variation,  the  result  Is  a 
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curve  wnich  falls  rapidly  to  an  asymptotic  ber.avior 
approaching  6  =  0.162.  Using  yager’s  correlation  with 
von  Kansan' s  approach  gives  a  prediction  for  ©  which  agrees 
closely  with  the  present  method  for  small  Reynolds  number 
and  becomes  asymptotic  to  9  =  0.217’  for  higher  Reynolds 
number.  Apparently  this  difference  is  due  to  the  fact 
that  the  logarithmic  law-cf-t  :e-wall  properly  accounts  for 
Reynolds  number  variation,  while  the  1/7  power  law  does  not. 

It  should  be  pointed  ou*  that  all  calculations  of  © 
in  the  present  study  assumed  an  initial  \alue  of  ©  =  0.226. 
This  initial  value  was  chosen  so  as  to  generate  a  monatonic 
curve  for  @  as  a  function  of  r.  It  was  found  that  perturbing 
this  initial  condition  caused  the  calculation  to  rapidly 
converge  to  the  curve  shown  in  Figure  3~7- 

3.^.2  Johnston's  Swept,  Forward -Facing  Step 

Johnston  (1970)  experimentally  investigated  the  flow 
over  a  two  inch  high  rectangular  step,  swept  at  forty  five 
degrees  to  the  main  flow  direction.  While  an  example  of 
a  fully  three-dimensional  situation,  this  flow  Is  not 
strictly  a  boundary  layer  due  to  the  large  variation  of 
pressure  in  the  nomal  direction.  For  the  purpose  of  the 
present  analysis  the  ’’ ;  reestreaai”  v*»j.ocitv  dist»*ibution  was 
determined  by  fitting  a  polynomial  co  J'.Ur.pton’ s  pressure 
coefficient  data.  Three  different  such  distributions  were 
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considered,  one  corresponding  to  the  wall  pressure  coefficient, 
one  to  the  vertically  averaged  pressure  coefficient,  and 
one  assuming  a  constant  "freestream"  velocity.  As  expected, 
use  of  the  average  pressure  coefficient  gave  the  best 
results  but  none  of  the  calculations  agreed  well  with  the 
data. 

Bradshaw  (1971)  also  analyzed  Johnston's  data  using  a 
finite  difference  scheme.  His  calculations,  using  the  wall 
pressure  coefficient,  showed  similar  behavior  to  the  present 
results.  In  order  to  overcome  this  difficulty  Bradshaw 
determined  the  local  value  of  the  ’’ freestream"  pressure 
gradient  from  tie  two-dimensional  pressure  coefficient 
data  given  by  Johnston.  With  this  additional  complication 
Bradshaw's  final  results  agreed  fairly  well  with  the  data 
except  in  zne  region  close  to  separation.  In  this  context, 
separation  refers  to  the  condition  where  the  surface  flow 
in  this  'boundary  layer’’  has  become  parallel  to  the  step. 

The  present,  analysis  follows  directly  from  equations 
(3-36)  and  (3-39)-  Making  the  assumption  that  the  sweep 
is  infinite,  and  hence  that  all  variations  parallel  to  the 
s :ep  are  zero,  relates  derivatives  in  the  s  and  n  directions 
by 


Also  the  metric  coefficients  h^  and  h?  are  equal  to  unity. 
These  simplifications  result  in  two  ordinary  differential 
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equations  for  \ ^  and  9  as  functions  of  s*  =  s/L.  The 
characteristic  length  L  is  taken  to  be  the  distance  from 
the  step  to  the  initial  point  of  the  calculation,  along 
the  ''f reestream  line'1. 

The  results  of  this  analysis  are  shown  on  Figures 
3-9  and  3-10-  Also,  Bradshaw's  finite  difference 
calculation  is  displayed  for  comparison.  On  Figure  3-9> 
it  can  be  seen  that  use  of  the  vertically  averaged  pressure 
coefficient  gives  good  agreement  with  the  skin  friction 
up  to  about  s*  -  0.8.  The  calculation  using  the  wall 
pressure  coefficient  is  low  while  use  of  a  constant  free 
stream  velocity  gives  results  which  are  much  too  high. 
Bradshaw's  results  are  also  in  good  agreement  but  extend 
as  far  as  s*  =  O.65.  None  of  these  predictions  agree  with 
the  measured  separation,  w.v_ch  occurs  at  about  s*  -  0.93* 

Figure  3-10  shows  a  comparison  between  the  present 
theory  and  Johnston's  data  for  P.  Again  the  average 
pressure  gradient  case  gives  the  best  agreement  but 
predicts  separation  too  early.  Bradshaw's  calculation  is 
somewhat  better  but  also  shows  early  separation. 

These  results  should  not  be  construed  as  indicating 
some  intrinsic  deficiency  in  the  present  theory,  since  this 
flow  is  not  truely  a  boundary  layer  flow.  Actually,  in 
order  to  accurately  predict  the  swept  step  flow  field  one 
should  use  the  full  Navier- Stokes  equations  and  not  the 
boundary  layer  equations. 


9U 


nfejtWJ 
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Figure  3.10.  COMPARISON  OF  PRESENT  THEORY  FOR 
THREE  PRESSURE  DISTRIBUTION  ASSUMPTIONS  WITH 
THE  WALL  STREAMLINE  IANGENT  DATA  OF  JOHNSTON 
(1970)  FOR  FLOW  NEAR  A  SWEPT  STEP. 
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3-->.3  iat  a  of  Francis  and  Pie',  e 

Francis  and  Fierce  (196?)  meas;'-:  ec  P,  the  tangent  of 
the  angle  between  the  total  surface  shear  stre' s  vector  and 
the  shear  stress  in  the  freestreac  direction,  in  a  rectan¬ 
gular  channel.  The  10  inch  wide  channel  consisted  of  a  18 
inch  straigh*  inlet  section  followed  by  a  curved  section 
of  either  a  55  inch  centerline  radius  of  curvature  (rei'erred 
to  as  series  5)  or  a  25  inch  centerline  radius  of  curvature 
(referred  to  as  series  2)  and  a  «8  inch  straight  discharge 
section.  Most  of  the  measurements  were  taken  on  the  center- 
line  but  a  few  of  the  measurements  were  taken  2  inches 
either  side  of  the  centerline.  The  Reynolds  number  was 
1.0^+0.03x10^  per  foot  for  both  series  5  and  series  2. 

The  freest  ream  velocity  variation  was  about  0-3  percent  per 
foot,  while  the  crossflow  velocity  variation  was  from  30 
to  pu  percent  per  foot  in  the  radial  direction.  This 
experiment  p’-oviries  the  possibility  of  examining  effects 
produced  by  a  crossflow  pressure  gradient  without  the 
influence  of  a  strong  pressure  gradient  in  the  freestream 
flow  direction. 

Equations  (3-38)  and  (3-39)  are  solved  with  the 
following  simplifications: 
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where  x*,  y*,  anc  ~c*  are  respectively  the  streamwise  flow 

direction,  the  direction  normal  to  the  st Teamwise  flow 

direction,  and  the  radius  of  curvature,  all  non-dimensior.al- 

ized  with  respect  to  the  length  of  the  channel.  This 

results  in  two  coupled  partial  differential  equations  for 

and  o  ss  functions  of  x*  and  y*.  I*  is  assumed  that 

\  and  ®  are  known  along  some  initial  line  x*  -  no  side 
'1  c 

boundary  conditions  need  be  known.  The  partial  derivatives 

occurring  in  these  equations  may  be  numerically  approximated 

in  a  number  of  ways.  For  illustrative  purposes,  they  were 

calculated  by  forward  differences.  Backward  or  central 

differences  ccuic  have  been  used  *'ust  as  eerily.  Also, 

if  more  accuracy  is  desired  it  is  possible  to  keep 
'1 

correcting  — ■ y  and  *  with  their  average  values  iteratively 
The  results  of  these  caiuclations  for  the  center  stream 
line  are  compared  to  the  measurec  values  cf  »  and  tt  the 
skin  friction  as  computed  by  Francis  and  Pierce  usirg  the 
Ludwieg  and  Tillmann  (1950;  shear  law,  as  shown  in  Figures 
3-11,  3-lf ,  3-13,  end  3~1‘**  Agreement  with  tne  data  over 
about  the  first  half  of  both  channels  is  qui’.e  good,  but 
then  the  theory  and  the  data  diverge.  The  skin  friction  as 
computed  by  the  present  theory  is  over  than  that  computed 
by  the  Ludwieg  and  Tillmann  law,  and  the  theoretical  values 
of  P  are  considerably  higher  than  the  measured  values.  Cne 
•possible  reason  for  these  discrepancies  is  that  Mager’s 
profile  may  not  adequately  represent  Francis  and  Pierce’s 
data.  Shanebrook  and  Hatch  (1970)  suggested  a  more 
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Figure  3-12.  COMPARISON  OF  THEORY  WITH  WALL 
STREAMLINE  TANGENT  DATA  OF  FRANCIS  AND 
PIERCE  (1967),  SERIES  5. 
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Figure  3.13.  COMPARISON  OF  PRESENT  THEORY  WITH 
EXPERIMENTAL  SKIN  FRICTION  AS  ESTIMATED  IN 
A  CURVED  DUCT  BY  FRANCIS  &  PIERCE  (1967) 
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sophisticated  family  of  hodograph  models  for  crossflow 
profiles.  Their  method  extended  Eichelbrenner' s  (1966) 
polynomial  representation  by  forcing  consecutive  higher 
order  derivatives  to  be  zero  at  the  wall  and  at  the 
boundary  layer  edge.  This  family  is  called  the 
family  of  hodograph  models  where  p^  is  the  number  of 
consecutive  zero  derivatives  at  the  wall  beginning  with 
the  ith  derivative  and  qj  is  the  number  of  consecutive 
zero  derivatives  at  the  boundary  layer  edge  beginning 
with  the  jth  derivative.  Besides  containing  many  terms, 
Shanebrook  and  Hatch’s  hodograph  models  also  contain  the 
angle  0  that  is  used  in  Johnston's  (1957)  triangular 
model,  Figure  3-3-  By  the  use  of  the  raomentun  integral 
equations  and  an  entrainment  equation.  Shanebrook  and 
Hatch  compared  their  hodograph  model  with  Mager's  velocity 
correlation  for  Francis  and  Pierce's  data.  Shanebrook 
and  Hatch's  model  shows  better  agreement  with  the  data, 
hut  the  same  trends  as  the  present  theory,  as  seen  in 
Figures  3-12  and  2-14.  However,  Nash  and  Patel  (1972) 
showed  that  Mager's  correlation  was  in  excellent  agreement 
with  the  crossflow  velocity  data  for  Francis  and  Pierce's 
experiment. 

Therefore,  it  is  felt  that  the  disagreement  between 
the  present  theory  and  Francis  and  Pierce's  data  is  not  the 
fault  of  the  assumed  crossflow  velocity  correlation.  Rather, 
it  is  believed,  as  Nash  and  Patel  pointed  out,  that  the 
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boundary  layer  in  Francis  and  Pierce’s  curved  channel  was 
dominated  by  the  effects  of  comer  flow  at  the  Junction  of 
the  side  walls  and  the  floor  of  the  channel-  These  effects 
were  not  accounted  for  in  the  present  analysis. 

3.4.4  Data  of  Klinksiek  and  Pierce 

As  is  pointed  out  in  the  introduction  to  this  chapter, 
Mager’s  velocity  correlation  is  not  capable  of  represent¬ 
ing  the  so-called  S-shaped  crossflow  velocity  profiles. 

Figure  3 “2-  However,  it  is  believed  that  the  present  theory 
will  give  accurate  results  even  if  the  crossflow  velocity 
profile  is  only  grossly  approximated;  as  long  as  ft  is  well 
predicted,  and  apparently  Mager' s  correlation  predicts  a 
accurately-  KLlnksiek  and  Pierce  (1970)  conducted  an 
experiment  in  which  they  obtained  S-shaped  crossflow 
velocity  profiles.  Their  experiments  were  conducted  in  a 
doubly  curved  channel  of  rectangular  cross  section.  An 
initial  straight  section  of  76.7  indie;,  was  followed  by  a 
60  degree  bend  with  a  25  inch  centerline  radius  of  curvature. 
This  in  turn  was  followed  by  a  12  inch  straight  section  and 
a  60  degree  bend  in  the  opposite  direction  with  a  55  inch 
centerline  radius  of  curvature.  The  same  equations  (3-36) 
and  (3-39)  tha*  were  used  to  compare  with  the  data  of 
Francis  and  Pierce  apply  to  this  situation.  The  results  are 
shown  in  Figures  3-15  and  3-16.  The  theoretical  values  of 
ft  are  considerably  above  the  data,  presumably  again  because 
of  side  wall  influence.  The  S-shaped  profiles  were  observed 
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Figure  3.15.  COMPARISON  OF  PRESENT  THEORY  WITH 


ESTIMATED  SKIN  FRICTION  IN  THE  CURVED  DUCT 
EXPERIMENT  OF  KIINKSIEK  &  PIERCE  (1970) 


STRAIGHT  .  25  INCH 

SECTION  I  RADIUS 


RADIUS 


STR.  55  INCH  NEGATIVE 
SEC  I  RADIUS  OF  CURVATURE 


ARC  LENGTH  S  -  Inches 


Figure  3.16.  COMPARISON  OF  PRESENT  THEORY  WITH  WALL 
STREAMLINE  TANGENT  DATA  OF  KLINKSIEK  &  PIERCE  (1971) 


in  the  lafer  curved  portion  of  the  channel  with  the 
inch  centerline  radius  of  curvature  (starting  at  or.1; 
inches  in  Pigu-es  3—15  ana  3-16)*  Bv  this  point  in  the 
calculation  the  theory  considerably  diverges  from  the  data. 
Thus  no  conclusion  can  be  drawn  about  the  general  use  of 
Mager's  correlation  in  situations  where  bilateral  shewing 
of  tne  crossflow  profile  is  p  esent. 

3*5  Summary 

The  approach  to  calculating  three-dimensional  skin 
friction  presented  in  this  chapter  seems  to  offer  great 
P'-omise*  It  is  an  integral  method  which  utilizes  "law- 
of-the-wall  velocity  correlations  and  results  in  two, 
coupled,  partial  differential  equations  having  the  stream- 
wise  shin  friction  and  the  angle  between  the  total  surface 
shear  stress  vector  and  the  shear  stress  in  the  freestream 
direction  as  *he  aependent  variables*  This  analysis 
contains  significantly  less  empiricism  than  "classical' 
integral  methods  and  is  computationally  much  simpler  than 
currently  available  differential  approaches- 

An  area  requiring  more  study  is  the  oroblem  of 
specifying  a  suitable  crossflow  velocity  correlation. 
Mager’s  hodograph  model  was  used  here  as  it  would  seem  to 
offer  the  oest  current  compromise  between  accuracy  and 
complication.  However,  it  is  not  capable  of  predicting 
crossflow  profiles  with  reversed  shewing;  and  this  must 
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be  regarded  as  a  deficiency  for  any  generally  applicable 
theory.  It  is  felt  that  the  development  of  a  more 
accurate  crossflow  model  must  wait  upon  the  availability  of 
more  accurate  profile  data. 

The  approach  developed  here  compares  extremely  well 
with  the  rotating  disk  data  of  Cham  and  Head  and  with  other 
predictions  of  this  flow.  The  relatively  poor  performance 
of  the  present  method  when  compared  to  the  data  of  Johnston, 
Francis  and  Pierce,  anc  Klinksi.ek  and  Pierce  should  not 
be  taken  as  an  indication  of  any  serious  deficiency  in  the 
theory.  Johnston’s  swept  step  flow  was  not  truly  a 
boundary  layer  due  to  the  presence  of  a  large  vertical 
pressure  gradient  and  the  Francis/Pierce  and  Kllrxksiek/ 
Pierce  data  sets  were  strongly  influenced  by  channel  corner 
flows.  Hius,  serious  evaluation  of  the  present  method  must 
await  the  availability  of  a  clean  set  of  experimental  data 
for  a  three-dimensional  boundary  layer  flow.  Such  a  set  of 
data  should  include  direct  measurements  of  skin  friction  in 
both  the  freestreaa  and  crossflow  directions  as  well  ac 
accurate  measurements  of  the  pressure  gradient  along  the 
free st re  am  line  and  velocity  profile  data. 
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APPENDIX  A 


A.l.  The  coefficients  for  use  in  the  heat  transfer 
differential  equations.  Eons- (1-29, 30),  are  defined  as: 
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Hie  various  partial  derivatives  which  are  needed  in  the  above 
expressions  may  be  evaluated  as  follows: 
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he  following  are  curve-.lt  expressions  for  the  coefficients  H. 
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In  order  to  use  these  curve-fit  expressions,  one  first 
evaluates  a,  0,  and  y  from  their  definitions,  Eqnt . (1-32), 

J. 

and  then  evaluates  5  from  the  two-diafeusional  law-of- 
the-wall,  Eqn. (1-33) -  The  only  other  parameter  in  the 
formulas  for  Is  u*,  which  is  directly  related  to 
skin  friction  by  the  relation  u*  *  xC^/T^)1^2. 
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A. 2.  The  coefficients  G,  for  use  in  the  three-dimensional 

i 

theory,  Eqns. (3-38, 39),  are  defined  as  follows: 
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Using  the  law-cf-the-wall  and  the  Mager  hodograph,  these 
integrals  may  be  evaluated  numerically  as  functions  of  the 

■A- 

parameters  X^,  og,  0,  and  5  .  The  following  curve-fits 
are  suggested  for  the  numerical  results: 
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It  is  these  curve-fits  which  w°re  used  to  compare  the 
present  theorv,  Eqns  .>  (3-38, 39)  ,  with  the  three-dimensional 
flew  experiments  of  Figures  . 9  through  3.16. 
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